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NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


TECHNICAL  NOTE  D-964 

THE  USE  OF  HIGHER-ORDER  DIFFERENCE  METHODS  IN  BEAM 
VIBRATION  ANALYSIS* 


By  Donald  T.  Greenwood** 


SUMMARY 


Simple  and  higher -order  difference  methods  for  the  solution  for  the  natural 
frequencies  of  vibration  of  a  uniform  beam  are  compared.  The  same  basic  higher- 
order  method  is  used  throughout  for  the  interior  cells,  but  three  different 
methods  of  boundary-condition  representation  are  given. 

Tables  and  graphs  of  the  error  in  mode  frequencies,  as  compared  with  a 
continuous  beam,  are  given  for  the  various  methods  as  a  function  of  the  number 
of  cells.  It  is  concluded  that  higher-order  methods  improve  accuracy  for  a 
given  number  of  cells,  with  essentially  no  change  in  the  quantity  of  com¬ 
puting  equipment  required. 


*Hie  work  reported  here  was  performed  under  NASA  Contract  No.  NsG-63-60,  ad¬ 
ministered  through  The  University  of  Michigan's  Office  of  Research  Adminstra- 
tion,  Ann  Arbor. 

**Associate  Professor  of  Aeronautical  and  Astronautical  Engineering,  The  Un¬ 
iversity  of  Michigan. 
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I .  INTRODUCTION 


One  of  the  most  commonly  encountered  problems  of  structural  dynamics  is 
that  of  the  lateral  vibrations  of  beams.  With  the  advent  on  a  large  scale  of 
analog  and  digital  computers  of  greatly  improved  speed  and  computational  ef¬ 
ficiency,  it  is  worthwhile  to  take  another  look  at  the  methods  available  for 
performing  these  computations.  This  report  is  principally  concerned  with  im¬ 
provements  in  analog  computer  methods  of  analysis  of  the  beam-vibration  prob¬ 
lem.  In  particular,  it  is  concerned  with  methods  applicable  to  the  electronic 
differential  analyzer.  Nevertheless,  it  will  be  seen  that  the  matrix  formula¬ 
tion  of  the  problem  and  most  of  the  results  are  directly  applicable  to  digital 
computation  as  well . 

Until  quite  recently,  all  attempts  to  use  analog  computers  in  the  solu¬ 
tion  of  partial  differential  equations  by  difference  techniques  were  re¬ 
stricted  to  simple  methods  with  accuracies  of  second  order  in  the  cell  size 
at  best.  In  1956,  Fisher^-  suggested  the  use  of  higher-order  differences  in 
the  solution  of  these  problems  on  the  analog  computer  and  indicated  some  of 
the  advantages  to  be  expected.  However,  he  was  not  concerned  in  his  examples 
with  the  beam  equation  and  did  not  study  the  problem  of  boundary-condition 
representation  using  higher-order  methods. 

This  report  considers  in  detail  the  adequate  representation  of  boundary 
conditions  when  higher-order  methods  are  used,  and  presents  the  results  of 
rather  extensive  computations  which  enable  one  to  obtain  a  feel  for  advantages 
and  disadvantages  of  each  method.  The  analysis  in  the  main  body  of  the  re¬ 
port  is  concerned  with  uniform  beams.  The  mode -frequency  errors  for  a  tapered 
cantilever  beam  are  considered  in  Appendix  B  and  comparisons  are  made  with 
the  continuous  beam. 

No  analysis  of  mode-shape  errors  has  been  included.  However,  spot  checks 
have  shown  these  errors  to  be  comparable  in  magnitude  with  the  frequency  errors. 

The  author  wishes  to  acknowledge  the  contributions  and  helpful  sugges¬ 
tions  of  R.  M.  Howe  of  the  Department  of  Aeronautical  and  Astronautical  Engi¬ 
neering,  particularly  with  respect  to  the  section  on  errors.  Thanks  are  also 
due  C.  K.  Shah,  who  did  most  of  the  numerical  computation,  and  J.  W.  Thatcher, 
who  aided  in  the  digital  computer  programming. 
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II .  SYMBOLS 


A  mode -amplitude  matrix 

ao>ai,etc.  constant  coefficients 

B  - [T]"1  [K] 


Bi>B2>etc . 


b 


Y>o ,bi,etc  • 


C 


C  i  ,  C2  ,etc. 


D 


E 


F 

f 

G 

H 

h 

I 

i 

K 

k 

‘9 

L 


difference  matrices 

an  element  of  matrix  B 

constant  coefficients 

capacitance,  constraint  matrix 

constants 

difference  matrix 

Young's  modulus 

constraint  force 

lateral  force  per  unit  length 

constant 

constant 

cell  size 

moment  of  inertia  of  cross  section,  unit  matrix 
mode  number 

system  matrix  [m]  1  [k ] 
stiffness  matrix 
bending  stiffness  of  a  cell 
inductance,  total  length  of  beam 


V 


M  bending  moment 

m  mass,  mass  matrix 

N  total  number  of  cells 

n  station  designation,  number  of  degrees  of  freedom 
R  resistance 

S  S  -  0  is  the  characteristic  equation 

T  -[K][B]  Also  the  time -dependent  factor  in  the  beam  deflection 
t  time 

V  shear  force 
x  position  along  the  beam 

Y  deflection  amplitude 

y  lateral  deflection 

P  eigenvalue  for  a  continuous  beam  of  unit  length 
A  a  small  perturbation  of  the  designated  quantity 
e  frequency  error 

C  dimensionless  lateral  deflection 

H  the  square  of  a  mode  frequency 

9  beam  slope 

A  eigenvalue  for  beam  of  length  N 

|  dimensionless  position  along  beam 

p  mass  per  unit  length 

T  dimensionless  time 

t 


_pr  OWO  CJ 


eigenvalue  for  difference  beam  of  unit  length 
indices  referring  to  row  and  column,  respectively 
reference  value 

transpose  of  the  designated  matrix 
computed  in  the  modal  (y)  coordinate  system 
approximately  equal  to 

amplitude  of  a  sinusoidally  varying  quantity,  determinant  of  a  matrix 
row  matrix,  or  usual  parenthesis  notation 
rectangular  matrix,  or  usual  brackets  notation 


column  matrix,  or  usual  braces  notation 
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III .  CONTINUOUS  BEAMS 


A.  THE  BEAM  EQUATION 

We  will  restrict  ourselves  to  the  study  of  the  lateral  vibrations  of  a 
thin  beam  where  linearity  is  assumed.  Ihe  partial  differential  equation  des 
cribing  the  motion  is 

£  i1  §)  * p  &  -  f<x’t)  (i> 

where 

x  =  position  along  beam 
y  =  lateral  deflection 

t  =  time 

El  =  bending  stiffness 

p  =  mass  per  unit  length 
f  =  lateral  force  per  unit  length 

It  can  be  seen  that  Eq.  (l)  is  of  fourth  order  in  x  and  second  order  in 
t.  For  our  purposes  it  is  often  convenient  to  write  four  equations  of  first 
order  in  x. 


where 


0  = 

& 

dx 

(2) 

M  = 

El 

ox 

(5) 

V  = 

dM 

(4) 

dx 

n 

-  ^  +  f(x,t) 
ox 

(5) 

0  =  slope  of  beam 

M  =  bending  moments 
V  =  shear  force 


-^ONVO  t) 
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Ihls  set  of  partial  differential  equations  is  equivalent  to  Eq.  (1).  We  will 
see  that  the  analog  computer  approach  to  the  problem  is  based  upon  this  latter 
set  of  equations  rather  than  the  original  fourth-order  equation. 


B.  SOLUTIONS  OF  THE  BEAM  EQUATION 

1.  General  Solution. — In  general,  the  parameters  p  and  El  are  functions 
of  x.  Initially,  however,  we  will  assume  that  each  is  constant.  Furthermore, 
we  will  assume  that  units  are  chosen  such  that  the  coefficients  are  unity. 

Our  principal  interest  is  in  the  solution  of  the  homogeneous  equation,  so  let 
us  set  f(x,t)  =  0.  Then  Eq.  (1)  becomes 


+  &Z  =  o  (6) 

dx4  dt2 


It  may  be  seen  that  essentially  the  same  equation  can  be  obtained  by  trans¬ 
forming  Eq.  (l)  to  dimensionless  form  and  again  setting  f(x,t)  =  0.  Let 

'  ■  E 

X 

L 

t_  fET 

L2  '  P 

where  L  is  the  length  of  the  beam.  Then 

+  S  -  0  (7) 

df;4  St2 

O  2 

Equation  (6)  can  be  solved  by  the  method  of  separation  of  variables. 

The  solution  is  of  the  form 


y 


Y(x)  T(t) 


(8) 
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where 


Y(x)  *  Ci  cos  \/p  x  +  C2  sin  •J  p  x  +  C3  cosh  s/p  x  +  C4  sinh'/p  x  (9) 


and 


T(t)  »  G  cos  pt  +  H  sin  Pt 


(10) 


Hie  solution  is  valid  only  for  certain  discrete  values  of  the  parameter  p. 

Ihe  constants  Ci,  C2,  C3,  and  C4  are  evaluated  from  the  conditions  at  the  ends 
of  the  beam.  Ihe  constants  G  and  H  are  evaluated  from  initial  conditions. 

Let  us  now  consider  the  equation  for  the  eigenvalue  P  for  various  end 
conditions . 

2.  Beam  Built-In  at  Both  Ends. — Hie  boundary  conditions  for  a  beam  of 
unit  length  built-in 'at  both  ends  are 


y(0)  *  y(l)  *  0  (11) 

9(0)  =  9(1)  =  0  (12) 

From  Eqs.  (9),  (ll) ,  and  (12),  we  find  that 

Ci  +  C3  =  0  (13) 

C2  +  C4  =  0  (14) 


(cos  "/p  -  coBh  Vp)Ci  +  (sin>/p  -  sinh  ,/p)c2  *  0  (15 ) 

(-  sin  '/p  -  sinh  v^pJCi  +  (cos  ■/p  -  cosh'/p)C2  *  0  (16) 


giving  the  characteristic  equation 


-fs-ONVO  tt 


fr  ON  VO  CJ 
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1  -  cos  v/p  cosh  s/p  -  0 


(17) 


and  also 


C2  cos  v/g  -  cosh  Vp 
Ci  sinh  \/"p  -  sin  Vp 


(18) 


The  first  four  eigenvalues  are 


Pi  =  22.37329 

Pa  *  61.67282 

Pa  -  120.9034 
p*  =  199.8594 


The  value  p  ■  0  has  no  practical  significance  in  this  case  and  is  omitted. 

The  ratio  C2/Ci  is  nearly  equal  to  -1  for  all  modes,  the  approximation 
being  particularly  close  for  modes  higher  than  the  first. 


It  will  be  helpful  in  the  later  development  to  express  the  form  of  the 
deflection  curve  as  a  power  series  in  x.  Using  a  series  expansion  for  the 
trigonometric  and  hyperbolic  functions  in  Eq.  (9)>  we  obtain 


Y(x) 


61 


pPx10  + 

10  i 


or 


Y(x)  * 


(20) 


Using  Eqs.  (2),  (3), 


and  (4),  we  can  obtain  the  other  oscillation  magnitudes: 


I®  I 


2CiVp 


x 


-fid 

2: 


pix8  p3xa 


51 


6: 


+ 


+ 


(21) 
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3«  The  Free-Free  Beam. — For  the  case  of  a  free-free  beam,  the  boundary 
conditions  are 


M(0)  =  M(l)  =  0 
V(0)  =  M(l)  =  0 


(24) 

(25) 


Substituting  into  Eq.  (9),  we  obtain 


-Cj.  +  C3  =  0 

(26) 

-C2  +  C4  =0 

(27) 

cos  ■J&j 

(Cl  +  (sinh  vTp  -  sin  v/^Cs  =  0 

(28) 

sin  v/p 

lex  +  (cosh  Vp  -  cos  'JlijCz  =  0 

(29) 

giving,  as  before, 


Cg  cos  >/p  -  cosh  Vg 
Cl  sinh  v/"p  -  sin  '/p 


(50) 


The 


eigenvalue  equation  is  the  same  as  for  the  built-in  case,  namely, 


1  -  cos  '/"p  cosh  Vp  =  0 


(31) 


Q  On  VO 
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Hie  value  3=0  has  physical  significance  in  this  case  but  since  the  mode  in¬ 
volves  no  elastic  deformation,  it  is  omitted.  The  remaining  eigenvalues  are 
identical  with  those  for  the  built-in  case. 


The  amplitude  of  the  motion  is 
Y(x)  -  2C 


lU  ♦...)♦  2c  (V5  ,♦  4? ♦  4?  ♦ ... 

4'  8.'  5.'  9! 


k:  8: 

giving  the  approximate  amplitudes 
Y(x)  =  2Ci(l  -  ■/£  x  + 


© 


(i  _vi  X  + ^ 

\  4!  51  8!  91  J 

=  2C*/e(i  +  ±  +  i*L.£*  +  .m\ 

\  3:  4:  ti  si  y 

|m  1  =  2clf/^-^  +  ^-i*I  +  ..>) 

Vi  31  61  71  y 

|V  I  2  2Ci3^  ^/i  x  -  ^ 


(32) 

(33) 

(34) 

(35) 

(36) 


If  we  compare  the  solution  for  the  free-free  beam  with  that  for  a  beam 

built-in  at  both  ends,  we  note  some  important  similarities.  We  have  seen  that 

the  nonzero  eigenvalues  are  identic al.  Beyond  that,  the  roles  of  the  deflec¬ 
tion  y  and  the  bending  moment  M  are  interchanged  in  the  two  cases  as  may  be 

seen  by  comparing  Eqs.  (20)  and  (35)  and  also  Eqs.  (22)  and  (33)*  Similarly, 

the  roles  of  the  slope  6  and  the  shear  force  V  are  interchanged.  This  sim¬ 
ilarity  will  hold  even  for  the  finite  difference  solution  which  will  be  ob¬ 
tained  later.  Therefore  it  will  not  be  necessary  to  carry  out  separate  calcu¬ 
lations  for  built-in  and  free-free  beams. 

4.  The  Cantilever  Beam. — The  boundary  conditions  in  this  case  are 


y(0)  =  0 

0(0)  =  0 

M(l)  =  0 

V(l)  =  0 


(37) 

(38) 

(39) 

(40) 


Applying  these  boundary  conditions  to  the  general  solution  given  by  Eq.  (9), 
we  obtain 


Ci  +  C3  =  0 

C2  +  C4  =  0 


(41) 

(42) 


!•»  WWiHWlVW:  JS*it  «*«»; 
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(cos  Vp  +  cosh  </p) Ci  +  (sin  Vp  +  8inh  Vp^Ca  ■  0 
^inh  Vp  -  sin  Vf^ Ci  +  (cos  Vp  +  cosh  V^C2  -  0 
From  Eqs .  (43)  and  (44 ',  we  obtain  the  eigenvalue  equation 


(43) 

(44) 


1  +  cos  Vp  cosh  Vp  -  0 


(45) 


and  also 


C2  =  .  cos  Vp  +  cosh  Vp 
sin  Vp  ■*■  sinh  Vp 


(46) 


the  ratio  being  approximately  -1  for  all  modes  higher  than  the  firBt. 
The  first  five  eigenvalues  are 


Pi  *  3 .516015 

fe  »  22.05449 

03  =  61.69721 

04  -  120.9019 
05  *  199.8595 


Note  that,  except  for  1  «  1,  the  value  of  0^  in  this  case  in  approximately 
equal  to  0i-i  for  the  built-in  or  free-free  case. 

The  amplitudes  are  given  by  the  equations 


«  CTWO  -=t- 


•r  owo  o 
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Each  of  these  expressions  contains  the  same  powers  of  x  as  was  obtained  for 
the  corresponding  built-in  case  in  Eqs.  (20)  to  (23)* 

5.  The  Pinned  Beam. —The  boundary  conditions  for  a  pinned  beam  sure 


y(o) 

II 

H 

II 

0 

(51) 

M(o) 

=  M(l)  =  0 

(52) 

Applying  these  conditions  to  the  solution  given  in  Eq.  (9),  we  obtain 


Ci  =  C3  =  C4  =  0 


(53) 


and  therefore 


Y(x)  =  Cz  sin  x 


(54) 


where  the  eigenvalues  are 


2,2 


3i  =  i2*' 


(i  -  1,  2,  3,  ..•) 


(55) 


Ik 


The  power  series  form  of  the  solution  is 


Y(x) 


(56) 


Also 


I®  I  - 
|m  | 
|V| 


c2  M  +  ££  .  £!£  , 

\  2!  41  61 

(Tvix  +  &i-^  + 

\  5i  51 

(l  +  M  -  ££  + 

V  21  41  V 


-  C2p 


a 

=  C2P2 


(57) 

(58) 

(59) 


The  eigenvalues  in  this  case  are  not  close  to  those  calculated  in  the 
previous  cases.  In  fact,  for  the  higher  modes  (large  i),  the  eigenvalues  tend 
toward  positions  midway  between  those  for  the  built-in,  free-free,  or  cant¬ 
ilever  beams . 


on  vo  a 
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IV.  DIFFERENCES  OF  ORDER  h2 


Previously  we  have  obtained  solutions  for  the  partial  differential  equa¬ 
tions  describing  a  continuous  uniform  beam.  Unfortunately,  the  analog  com¬ 
puter  cannot  solve  partial  differential  equations  directly  but  must  solve  the 
set  of  ordinary  differential  equations  which  are  obtained  by  finite  difference 
approximations  to  the  spatial  derivative.  Thus  we  will  obtain  solution  ampli¬ 
tudes  at  a  discrete  set  of  points  rather  than  continuously  as  a  function  of  x. 
At  the  given  points,  however,  the  finite  difference  solution  should  closely 
approximate  the  solution  for  the  continuous  beam. 


A.  FINITE  DIFFERENCE  EQUATIONS 


An  analog  computer  approach  to  the  beam-vibration  problem  using  finite 
difference  equations  of  order  ha  has  been  given  by  Howe  and  Howe**  and  others 
The  procedure  in  each  case  is  based  essentially  upon  Eq.  (60). 


■yn  +  yn+i  h2  d3y 
h  2k  Ux3 


n+^ 


(60) 


The  interval  between  stations  is  designated  by  h.  Subscripts  refer  to  the  sta¬ 
tion  at  which  y  (or  one  of  its  derivatives)  is  calculated.  The  last  term  is 
the  first  error  term  and  is  not  included  in  the  computer  mechanization.  It  is 
given  to  provide  an  estimate  of  the  accuracy  of  the  approximation.  Equation 
(60)  can  be  derived  by  writing  the  Taylor  expansions  about  station  n-H^  for  yn 
and  yn+1  and  solving  for  after  eliminating 

°x | n+5 

Using  the  approximation  of  Eq.  (60)  and  again  assuming  that  El  and  p  are 
unity,  the  beam  equations  (2),  (5),  (4),  and  (5)  can  be  written  in  difference 
form: 
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(65) 
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Note  that  displacements  and  bending  moments  are  calculated  at  Integer 
stations  whereas  the  slopes  and  shear  forces  are  calculated  at  half -Integer 
stations.  Zf  one  calculates  all  quantities  only  at  Integer  stations,  the  re¬ 
sult  is  to  double  the  effective  cell  length  and  thereby  to  quadruple  the  first 
error  term.  Therefore  we  will  use  "staggered"  stations  when  finite  difference 
approximations  to  first  derivatives  are  used. 

The  analog  computer  circuit  representing  a  typical  internal  cell  is 
shown  in  Fig.  1.  It  requires  five  amplifiers  per  cell  if  one  reverses  signs 
at  adjacent  cells. 

Another  approach  to  the  problem  is  to  approximate  tne  second  spatial 
derivative  rather  than  the  first  derivative.  In  this  case  the  finite  dif¬ 
ference  approximation  is 


& 

dx2 
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yn-i  -  2yn  +  yn+i  hs 

h2  “  12  dx4 


(65) 


and  the  corresponding  beam  equations  solved  on  the  computer  are 


Mn 


yn-i  -  gyn  +  yn+i 

h2 


(66) 


■Mn-i  +  2Mn  -  Mn+1 


(67) 


The  analog  computer  circuit  for  this  case  requires  only  3  amplifiers  per 
dell  and  is  shown  in  Fig.  2.  This  circuit  is  subject  to  exactly  the  same 
finite  difference  errors  as  the  5 -amplifier  circuit  and,  in  fact,  can  be  ob¬ 
tained  as  a  direct  reduction  of  that  circuit. 


B.  PHYSICAL  INTERPRETATION 

There  are  two  ways  in  which  one  can  give  physical  meaning  to  finite  dif¬ 
ference  approximations  to  the  beam  equations. 
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The  first  method  is  to  consider  a  finite  difference  expression  as  an 
approximation  to  a  spatial  derivative  in  the  continuous  beam  at  a  given  point. 
For  example,  the  right-hand  side  of  Eq.  (64)  approximates  -  2jU  ,  which,  in 
accordance  with  Eq.  (5)  for  a  continuous  beam,  is  proportional  to  the  accelera¬ 
tion  yn  of  a  local  differential  element. 

The  second  method  is  to  think  of  the  set  of  finite  difference  equations 
as  representing  exactly  a  lumped  physical  system  whose  response  approximates 
that  of  the  continuous  beam.  The  calculated  values  of  accelerations,  shear 
forces,  etc.,  at  the  various  stations  have  their  real  counterparts  in  the 
lumped  physical  system.  This  approach  is  helpful  in  gaining  physical  insight 
into  the  meaning  of  the  equations,  particularly  for  simple  differences.  Two 
lumped  physical  representations  of  this  finite  difference  approximation  to 
the  beam  equations  are  shown  in  Fig.  3* 

In  Fig.  3&  "the  beam  is  approximated  by  a  series  of  massless  levers  con¬ 
nected  by  pin  Joints.  At  each  Joint  a  spring  produces  a  moment  (corresponding 
to  the  bending  moment)  which  is  proportional  to  the  difference  in  slope  angle 
9  of  the  adjacent  levers.  The  mass  of  the  beam  is  lumped  into  point  masses 
at  the  Joints.  The  deflection  y  is  measured  at  the  Joints  while  the  slope  9 
is  measured  at  the  midpoint  of  each  lever. 

Figure  3b  shows  an  electrical  circuit  analog  of  the  lumped  mechanical 
system  of  Fig.  3a-  It  also  Is  described  by  finite  difference  Eqs.  (6l)  to 
(64).  The  displacement  velocities  are  represented  by  voltages  and  the  shear 
forces  and  bending  moments  are  represented  by  currents.  The  lumped  bending 
stiffnesses  are  represented  by  inductors,  while  the  lumped  masses  are  repre¬ 
sented  by  capacitors.  The  transformers,  each  of  turns  ratio  1  to  h,  perform 
the  coordinate  transformations  relating  slopes  to  deflections,  and  can  be 
considered  as  the  analog  of  the  massless  levers. 


C.  BOUNDARY  CONDITIONS 

We  have  previously  given  the  equations  stating  the  boundary  conditions 
on  a  continuous  beam  for  the  various  cases  under  consideration.  We  must  now 
write  corresponding  equations  for  the  finite  difference  beam.  First,  however, 
we  must  decide  whether  to  end  the  beam  at  a  point  where  lateral  deflection 
and  bending  moment  are  computed  (the  "integer"  stations),  or  at  a  point  where 
the  slope  and  shear  force  are  computed  (the  "half -integer"  stations). 

Experience  has  shown^  that  higher  accuracy  can  be  obtained  for  a  given 
amount  of  computation  (or  degrees  of  freedom)  if  the  end  occurs  ht  a  half- 
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integer  station  for  built-in  or  free  ends.  For  the  pinned  case  it  is  best  to 
end  the  beam  at  an  integer  station. 

1.  Built-In  End. —Previously  we  have  seen  that  the  displacement  y  and 
slope  8  are  zero  at  a  built-in  end  of  a  beam.  Suppose,  for  example,  that  the 
beam  ends  at  station  n+£.  We  must  approximate  the  deflection  at  a  half¬ 
integer  station  where  it  is  not  usually  defined.  Using  the  Taylor  expansion 
approach  and  keeping  only  the  first  error  term,  we  find  that 

*«*  ‘  Ky»  *  y»«)  '  r  & 


Therefore,  since  yn+£  and  8n+^  are  zero,  we  see  from  Eqs.  (6l)  and  (68)  that 

yn  =  yn+1  =  o  (69) 

implying  that  there  is  no  lateral  displacement  at  the  integer  station  adjacent 
to  a  built-in  end. 

2.  Free  End — At  a  free  end  we  find  that  the  bending  moment  M  and  shear 
force  V  are  zero.  Using  reasoning  similar  to  that  of  the  previous  case,  we 
find  that  for  a  beam  ending  at  station  n+£ 


«n  =  Ku«  "  0 


(70) 


3-  Pinned  End . —For  a  pinned  end  at  station  n  the  boundary  conditions 
are  quite  straightforward,  namely, 


II 

o 

(71) 

Mn  -  0 

(72) 

In  case  the  pinned  end  occurs  at  station  n+£,  the  boundary  conditions 
are 
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-yn+1 

(73) 

Mn  * 

-Mn+i 

(74) 

4.  The  Use  of  Symmetry.— For  vibration  problems  of  a  uniform  beam  where 
the  same  boundary  condition  applies  at  each  end,  one  can  make  use  of  symmetry 
to  reduce  the  required  amount  of  computation.  All  the  natural  modes  can  be 
classified  according  to  whether  the  deflection  curve  has  even  or  odd  symmetry 
about  the  midpoint  of  the  beam.  Consider,  for  example,  a  beam  that  is  built- 
in  at  both  ends.  If  one  arranges  the  modes  in  order  of  increasing  frequency, 
then  modes  1,  3,  5 ,  etc.,  will  have  mode  shapes  with  even  symmetry  about  the 
midpoint,  whereas  modes  2,  4,  6,  etc.,  will  exhibit  odd  symmetry  about  the 
midpoint.  (See  Fig.  4.)  In  either  event,  however,  one  need  analyze  only  one- 
half  of  the  beam  since  the  motion  of  the  other  half  can  be  deduced  from  sym¬ 
metry. 

For  the  case  of  even  symmetry  of  the  deflection  curve,  the  bending  mo¬ 
ment  curve  will  also  show  even  symmetry,  whereas  the  slope  and  shear  force 
curves  will  have  odd  symmetry.  Conversely,  for  odd  symmetry  in  y  and  M, there 
will  be  even  symmetry  in  9  and  V. 

The  above  symmetry  assumptions  sure  exact  and  therefore  the  same  fre¬ 
quencies  and  mode  shapes  result  as  for  the  case  where  the  equations  are  writ¬ 
ten  for  the  complete  beam. 


D.  ANALOG  COMPUTER  CIRCUITS 

1.  Cantilever  Beam. — The  analog  computer  circuit  for  the  cantilever  beam 
is  useful  in  illustrating  built-in  and  free  end  conditions .  Two  circuits  are 
shown  in  Fig.  5-  The  first  circuit  requires  3  amplifiers  per  cell,  the  sec¬ 
ond,  5  amplifiers  per  cell.  These  are  nominal  values,  however,  and  amplifiers 
can  be  saved  in  representing  the  end  cells  because  of  the  boundary  conditions. 
Thus  it  can  be  seen  that  the  actual  numbers  of  amplifiers  required  to  repre¬ 
sent  a  four-cell  cantilever  beam  are  9  and  15,  respectively.  Cells  can  be 
added  by  adding  standard  circuits  such  as  those  in  Fig.  1  or  Fig.  2,  repre¬ 
senting  interior  cells. 

2.  Built-In  Beam.— The  computer  circuits  for  a  six-cell  beam  that  is 
built-in  at  both  ends  are  shown  in  Fig.  6.  Note  that  only  one-half  of  the 
beam  is  actually  represented  on  the  computer;  the  motion  of  the  other  half  is 
inferred  from  symmetry.  Circuits  are  shown  for  even  or  odd  symmetry  about  the 
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center  of  the  span.  In  each  case  the  center  occurs  at  a  half-integer  station. 
Circuits  for  the  case  where  the  center  occurs  at  an  integer  station  could 
have  been  derived  hy  a  similar  process.  This  situation  would  arise  for  a 
built-in  beam  with  an  odd  number  of  cells ,  assuming,  of  course,  that  the  ends 
occur  at  half -integer  stations. 

3«  Pinned  Beam. —The  computer  circuits  representing  a  six-cell  beam 
with  both  ends  pinned  are  shown  in  Fig.  7*  The  circuit  for  the  case  of  odd 
symmetry  about  the  center.  Fig.  7(b),  is  particularly  simple  because  of  the 
additional  symmetry  about  the  quarter  points  along  the  span. 

The  circuits  to  be  used  in  the  higher-order  methods  that  we  consider  are 
quite  similar  in  general  form  to  those  of  FigB.  5(b),  6,  and  7*  However, 
most  amplifiers  will  have  four  rather  them  two  inputs. 


E.  MATRIX  FORMULATION  OF  THE  DIFFERENCE  EQUATIONS 

In  the  emalysis  of  specific  cases  of  finite  difference  beams,  it  is  con¬ 
venient  to  write  the  equations  in  matrix  form.  It  is  particularly  important 
to  use  this  formulation  if  digital  computations  are  to  be  performed. 

In  this  report  the  beam  equations  sure,  in  general,  written  as  four  first- 
order  difference  equations.  Using  matrix  notation,  they  are  as  follows: 
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Although  the  B  matrices  are  rectangular  in  general,  the  K  matrix  must  be  square. 
The  natural  frequencies  are  obtained  from  the  determinants!  equation 


|K  -  X2l|  -  0 


(81) 


where  X  is  the  natural  frequency  in  rad/sec  and  I  is  the  unit  matrix.  For  the 
case  where  the  K  matrix  is  n  x  n,  there  will  be  n  values  of  X2  (roots)  cor¬ 
responding  to  the  n  natural  frequencies  of  the  system. 

As  an  example,  consider  the  case  of  a  six-cell  beam  that  is  built-in  at 
both  ends.  We  will  write  the  equations  for  modes  with  even  symmetry  at  the 
center.  As  before,  El  and  p  are  assumed  to  be  unity,  ftie  matrix  equations 
are 
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Assuming  a  unit  cell  size  (h 
tions,  we  obtain 


1)  and  performing  the  matrix  multiplica- 
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(86) 


According  to  Eq.  (8l),  the  roots  are  obtained  from 


(6-A2)  -3 
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giving  the  characteristic  equation 


\4  -  8a2  +  3  =  0 


(88) 


The  roots  are 


*!,a  =  0.}9khk9,  7*60555 
or 

Ai,2  =  0.628052,  2.75782 


These  are  the  actual  frequencies  for  the  modes  with  even  symmetry  and 
h  =  1.  For  our  purposes,  however,  it  is  convenient  to  assume  that  the  toted 
length  of  the  beam  is  unity.  This  implies  that  h  =  l/N,  where  N  is  the  number 
of  cells,  and  results  in  multiplying  each  frequency  by  N2.  Denoting  this 
normalized  frequency  by  00,  we  obtain 

u>i  =  22.6099 

a>s  =  99.2815 
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The  same  basic  matrix  formulation  will  be  used  in  the  computations  through¬ 
out  the  report.  Hie  use  of  higher-order  differences  will  add  to  the  numerical 
complexity  but  will  not  change  the  form  of  the  equations. 


F.  RESULTS 

The  results  of  computations  of  the  normalized  frequency  parameter  cn  using 
differences  of  order  h2  are  summarized  in  Tables  II  and  III.  The  computations 
are  for  built-in  (or  free-free)  and  cantilever  beams  for  various  numbers  of 
cells.  The  finite  difference  error  is  plotted  against  the  cell  number  N  in 
Figs.  13  and  14.  These  results  are  essentially  the  same  as  those  given  pre¬ 
viously  by  Howe  and  Howe^  and  Michie,7  except  that  the  higher  mode  fre¬ 
quencies  are  also  included. 


G.  ORTHOGONALITY 

We  have  seen  previously  that  the  equations  of  motion  can  be  put  in  the 

form 


(y)  +  [K]  (y)  =  0 


(79) 


We  can  also  write 


[m]  (y)  +  [k]  (y)  =  0 


(89) 


where  m  and  k  are  the  mass  and  stiffness  matrices,  respectively.  Comparing 
Eqs.  (79)  and  (89) ,  we  see  that 


[K]  =  [m]  [k] 


(90) 


Now,  for  the  case  of  the  first-order  difference  equations  under  considera¬ 
tion  here,  both  the  m  and  k  matrices  are  symmetric.  In  fact,  the  m  matrix  is 
diagonal.  It  is  always  possible  in  these  linear  bilateral  undamped  systems 
(whose  kinetic  and  potential  energies  can  be  written  as  positive -definite 
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quadratic  forms)  to  find  a  coordinate  transformation  which  will  simultane¬ 
ously  diagonalize  the  m  and  k  matrices.  In  this  case  the  natural  modes  of  the 
system  are  said  to  be  orthogonal. 

Ihe  orthogonality  of  the  modes  is  not  surprising  since  we  have  already- 
seen  (Fig.  3)  that  these  equations  represent  exactly  a  lumped  mechanical  or 
electrical  system  composed  of  linear  bilateral  elements.  However,  when  we 
discuss  higher-order  difference  methods,  we  will  find  that  the  k  matrix  is 
not  always  symmetric,  in  which  case  the  modes  are  not  orthogonal. 
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v.  DIFFERENCES  OF  ORDER  h4 


We  have  given  a  brief  review  of  the  calculation  of  eigenvalues  for  a  uni¬ 
form  beam,  using  simple  finite  difference  approximations  to  the  spatial  deriv¬ 
atives.  In  this  way  we  have  obtained  some  physical  feeling  for  the  main  fea¬ 
tures  of  this  method  and  perhaps  some  insight  into  the  more  promising  approaches 
to  the  use  of  higher-order  differences.  In  the  remainder  of  the  report  we 
will  consider  several  higher -order  methods  and  indicate  some  of  the  advantages 
as  well  as  problems  in  their  use. 


A.  CHOICE  OF  THE  BASIC  DIFFERENCE  EQUATION 

Even  before  one  gets  into  the  questions  concerning  the  representation  of 
boundary  conditions,  there  are  a  number  of  possible  choices  one  can  make  con¬ 
cerning  the  basic  higher-order  difference  equation  to  be  used.  First  is  the 
question  whether  the  approximation  ought  to  be  with  respect  to  the  first,  sec¬ 
ond,  or  possibly  the  fourth  spatial  derivative.  We  can  immediately  reject 
the  latter  possibility  on  the  grounds  that  it  is  quite  sensitive  to  component 
errors  in  the  analog  com outer  circuit.® >9  Also  there  is  the  disadvantage  that 
other  variables  of  interest  such  as  ©,  M,  and  V  are  not  directly  available 
even  though  three  amplifiers  are  required  per  cell. 

The  choice  between  approximating  two  second-order  equations  or  four  first- 
order  equations  is  more  difficult.  Let  us  consider  each  possibility  in  turn. 

The  most  obvious  approximation  of  order  h4  for  the  second  derivative  with- 
respect  to  x  is  given  by  the  central  difference  equation 


~yh-g  +  ^®yn-i  ~  3Qyn  *  ^-®yn+i  ~  yn+2 
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Using  this  approximation,  the  equations  set  up  on  the  computer  are  of  the  form 
Mjj  *  ryn-2  +  l®Fn-i  "  ^Oyn  +  ^Fn+i  “  (92) 
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where  unity  values  are  assumed  for  El  and  p.  Ihe  analog  computer  circuit  re¬ 
quires  three  amplifiers  per  cell  and  up  to  five  inputs  per  amplifier.  (See 
Fig.  8.) 

Consider  now  some  of  the  higher-order  approximations  to 
central  difference  approximation  of  order  h4  is 
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Ihe  resulting  beam  equations  are 
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where  El  and  p  are  assumed  to  be  of  unit  value.  The  analog  computer  circuit 
requires  nine  amplifiers  per  cell  in  this  case  because  each  of  the  four  output 
variables  must  be  generated  with  both  signs.  (Each  cell  would  require  2 
integrators, >3  summers,  and  4  inverters.) 

Another  possibility  is  to  try  a  higher-order  approximation  using  stag¬ 
gered  stations,  i.e.,  calculate  the  derivatives  at  the  midpoints  between  sta¬ 
tions  where  the  function  is  defined.  A  central  difference  approximation  of 
this  sort  is  the  following: 
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This  approximation  is  seen  to  have  a  smaller  error  term  than  that  for 
Eq.  (9*0*  Furthermore,  it  will  be  seen  that  less  analog  computer  equipment 
is  required  to  implement  it. 

Using  the  approximation  of  Eq.  (99) >  the  beam  equations  are 

®n+£  =  "  27yn  +  2Tyn+i  "  yn+2^ 

Mj,  -  —  ( 6  a-276  i  +  276  ,  -  6  ^ 

24h  v-i  n-^  n+^  n+y 

Vj  -  ' 21M"  *  27M"*‘ ' "»«) 

y»  ■  sk  (-M  *  27v-i  ■  27v* '  v»*t) 

where,  as  before,  unit  values  are  assumed  for  El  and  p.  Note  that  y  and  M  are 
calculated  at  integer  stations  whereas  6  and  V  are  calculated  at  half -integer 
stations . 

The  analog  computer  circuit  representation  of  Eqs.  (100)  to  (103)  re¬ 
quires  5  amplifiers  per  cell  with  up  to  4  inputs  per  amplifier.  The  circuit 
for  a  typical  cell  is  shown  in  Fig.  9.  The  signs  of  amplifier  outputs  alter¬ 
nate  on  successive  cells. 

In  spite  of  the  fact  that  the  circuit  of  Fig.  9  requires  5  amplifiers  per 
cell  compared  to  3  amplifiers  per  cell  for  the  circuit  of  Fig.  8,  the  computa¬ 
tions  in  this  report  are  based  on  the  beam  equations  as  given  by  Eqs.  (100) 
to  (103).  The  reasons  for  this  choice  are  (1)  the  greater  accuracy  of  the 
basic  difference  equation,  (2)  greater  ease  and  flexibility  in  applying  bound¬ 
ary  conditions  and  forcing  functions,  and  (3)  less  sensitivity  to  computer 
component  errorB. 


B.  BOUNDARY  CONDITIONS  AND  RESULTS 

Earlier  we  saw  that  the  representation  of  boundary  conditions  was  relatively 
straightforward  for  the  case  of  difference  equations  of  order  h2.  When  one 
uses  higher-order  differences,  the  problem  of  applying  boundary  conditions  be- 
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comes  more  complex.  Many  approaches  are  possible  and  the  accuracy  of  the  re¬ 
sults  Is  strongly  influenced  by  the  method  chosen.  Much  of  the  remainder  of 
this  report  will  be  concerned  with  a  study  of  the  most  promising  methods  of 
applying  boundary  conditions  when  using  higher-order  differences  in  beam- 
vibration  problems. 

It  can  be  seen  that  if  one  uses  the  basic  central  difference  approxima¬ 
tion  of  Eq.  (99) »  difficulties  arise  as  one  approaches  the  end  of  the  beam  be¬ 
cause  the  needed  stations  beyond  the  end  of  the  beam  do  not  exist  physically. 
Consequently,  one  must  make  assumptions  concerning  virtual  stations  or  else 
one  must  use  unbalanced  or  one-sided  formulas. 

1.  Pinned  Beam. — An  end  condition  which  is  easy  to  implement  is  the 
pinned  end.  In  this  case  one  can  extend  the  beam  by  assuming  symmetry  and  use 
these  virtual  stations  in  the  standard  central  difference  formula  that  we  have 
chosen.  Odd  symmetry  is  assumed  for  the  displacement  y  and  bending  moment  M, 
whereas  even  symmetry  is  assumed  for  the  slope  ©  and  the  shear  force  V.  These 
assumptions  are  entirely  in  keeping  with  the  actual  facts  for  the  continuous 
case  and  no  error  is  introduced  by  this  means.  Thus  any  error  in  the  fre¬ 
quency  calculations  is  due  to  the  basic  finite  difference  approximation  rather 
than  the  method  of  applying  boundary  conditions. 

The  results  of  frequency  error  calculations  for  a  beam  pinned  at  both 
ends  are  shown  in  Fig.  21.  These  curves  also  represent  the  error  inherent  in 
the  basic  central  difference  approximation.  They  indicate  that  the  calculated 
frequency  tends  to  be  low  with  the  error  decreasing  roughly  as  h4  for  in¬ 
creasing  numbers  of  cells. 

It  might  be  well  to  note  that  there  is  another  end  condition  which  can 
be  represented  without  error.  This  is  the  case  for  which  y  and  M  have  even 
symmetry  and  ©  and  V  have  odd  symmetry  about  the  end. 

2 .  Built-In  Beam. — The  two  cases  mentioned  above  are  the  only  ones  for 
which  we  can  use  the  basic  central  difference  equation  in  a  straightforward 
manner  for  virtual  stations  beyond  the  end  of  the  physical  beam.  All  the  other 
boundary  conditions  require  a  different  approach.  One  possibility  is  to  use 
one-sided  approximations  of  order  h4 .  This  means  that  the  input  data  for  cal¬ 
culating  the  spatial  derivative  at  a  point  do  not  come  from  stations  placed 
symmetrically  with  respect  to  that  point.  Rather,  a  majority  or  perhaps  all 

of  the  stations  sure  located  on  one  side  of  the  point  at  which  the  derivative 
is  calculated. 
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As  was  the  case  when  we  used  simple  differences  (error  of  order  h2)  we 
will  again  end  the  team  at  a  half -Integer  station.  Dms  9  and  V  are  normally 
calculated  at  the  end,  whereas  y  and  M  are  calculated  at  the  Integer  station 
that  Is  located  h/2  Inboard.  However,  since  the  boundary  conditions  may 
Involve  the  values  of  y  or  M  at  the  end  of  the  beam,  we  must  develop  suit¬ 
able  approximations  to  their  values  at  these  points. 

In  this  section  we  will  use  difference  approximations  of  order  h4,  re¬ 
quiring  input  data  only  from  stations  on  the  actual  beam.  No  symmetry 
assumptions  are  made  at  the  ends.  The  approximations  actually  used  in  the 
development  are  given  in  Table  I. 

Now  let  us  consider  the  case  of  a  beam  that  is  built-in  at  station 
l/2  located  at  its  left  end.  Using  one-sided  or  central  difference  approx¬ 
imations  as  required,  we  obtain 
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But  from  Table  I  we  note  that 
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where  the  error  terms  of  order  h4  and  h5  are  zero  for  the  case  of  a  uniform 
beam.  This  can  be  seen'  from  Eq.  (l)  or  (6)  in  conjunction  with  the  boundary 
conditions  yA  »  0  and  0^  =  0.  Also,  it  follows  directly  from  the  series  form 
2  2 

of  the  solution  as  given  by  Eq.  (19). 


Omitting  terms  of  order  h4  or  higher,  we  find  that 
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Thus  we  find  that  we  can  omit  yi  as  a  coordinate  in  our  analysis  because 
of  the  constraint  expressed  in  Eq.  (105)*  This  was  also  true  when  simple  ap¬ 
proximations  of  order  h2  were  used,  but  in  that  case  yx  was  equal  to  zero. 

In  a  similar  manner,  using  the  approximations  of  Table  I  and  Eq.  (99) , 
we  can  write  the  remaining  beam  equations  for  the  case  of  a  built-in  end. 
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Let  us  apply  Eqs.  (106)  to  (109)  for  the  specific  case  of  a  6-cell  beam 
that  is  built-in  at  both  ends.  Only  the  modes  with  even  symmetry  about  the 
center  will  be  considered.  Ihe  equations  of  motion  are 
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Letting  h  =  1  and  applying  Eqs.  (80)  and  (8l),  we  obtain 


(5.22099-*2)  -2.77555 
-4.29587  (2.85815.x2) 
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giving  the  characteristic  equation 


-  8.0791^X2  +  3.00462  -  0 


(115) 


and  roots 


Xi>2  *  0.625142,  2.77278 


which  can  be  normalized  to 


OJi  *  22.5051 
ere  =  99.8203 


Comparing  these  results  with  those  for  the  case  where  simple  differences 
were  used,  one  finds  that  the  h4  approximation  gives  somewhat  better  accuracy 
for  the  first  mode  and  approximately  the  same  results  for  the  higher  mode. 
This  improvement  due  to  using  an  h4  approximation  seems  to  be  rather  slight. 
However,  if  one  looks  at  the  accuracy  of  a  given  mode  as  the  number  of  cells 
is  increased,  the  h4  approximation  converges  toward  the  correct  result  more 
rapidly,  as  it  must,  of  course,  in  the  limit.  See  Table  II  and  Fig.  15  for 
a  summary  of  the  results . 


3.  Free  Beam. — Previously  we  have  seen  that  the  representation  of  a  free 
end  is  obtained  from  that  for  a  built-in  end  by  interchanging  the  roles  of  y 
and  M  and  also  those  of  0  and  V.  Furthermore,  we  found,  for  a  continuous  beam 
that  the  nonzero  mode  frequencies  are  identical  for  the  cases  where  both  ends 
are  built-in  or  both  ends  are  free.  This  is  also  true  when  finite  difference 
methods  are  used.  To  illustrate  this,  let  us  calculate  the  mode  frequencies 
(even  symmetry  only)  of  a  6-cell  beam  that  is  free  at  both  ends. 


The  equations  of  motion  are 


-93  207  -114 

f 

yi 

2 

9a 

2 

„  =  — 

>  24h 

-22  18  4 

< 

y2  - 

% 

1  -27  26 

y3 

2 

-»  “ 

—  + 

h  4 

t 


(116) 


■P"  o\\o  t> 


P~  OWO  tJ 


55 


Me 

> 

M3 
>  < 


-27 


1 


f  • 

27 

®1 

2 

®a 

-26 

p. 

®5 . 

2 

k*  / 

w 

2! 

>  *  1 

4800  0 

' 

f  N 
M2 

>  Z> 

5400h 

-6025  5841 

M3 

t  . 

✓  > 
yi 

-17  -8 " 

f  \ 

va 

y s 

(  _  i- 

27  -27 

1 

vf 

y3 

k  > 

24h 

-1  26 

*5. 

2 

*  4 

(117) 


(118) 


(119) 


Performing  the  matrix  multiplications  for  the  case  where  h  =  1,  we  obtain 
from  Eq.  (8l)  that 


(0.20619-X2) 

-2.16868 

1.23959 


-0.00031 
( 5. 46771 -X2) 
-3.64483 


-0.20588 

-3.29903 

(2. 40524 -X2) 


0 


(120) 


giving  the  characteristic  equation 


X2(X4  -  8. 07914 X2  +  3.00462)  =  0 


(121) 


Comparing  Eqs.  (115)  and.  (121),  we  see  that  the  nonzero  roots  are  identical 
even  though  the  K  matrices  are  quite  different. 

All  computations  in  this  report  applicable  to  either  a  built-in  or  free 
beam  were  actually  set  up  considering  the  built-in  case,  using  equations 
similar  to  (110)  to  (113). 

4.  Cantilever  Beam. —One  can  use  the  methods  developed  for  obtaining 
h4  accuracy  in  the  representation  of  a  built-in  end  or  a  free  end  to  write  the 
equations  of  motion  for  a  cantilever  beam.  The  minimum  number  of  cells  re¬ 
quired  for  h4  accuracy  is  five.  Other  methods  to  be  considered  later  will  be 


5 


applicable  with  fewer  cells. 

As  an  example  of  a  typical  calculation,  consider  a  six-cell  cantilever 
beam.  Hie  equations  of  motion  are 
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Hie  large  numbers  in  Eqs.  (122)  and  (124)  result  directly  from  using  the 
one-sided  approximations  of  order  h4  as  given  in  part  A  of  Table  I. 

Note  the  similarities  in  the  matrices  of  Eqs.  (122)  and  (124).  Reading 
from  left  to  right,  beginning  with  the  first  row  of  (122),  one  sees  the  neg¬ 
ative  of  the  numbers  obtained  in  reading  from  right  to  left  beginning  with 
the  last  row  of  (124).  A  similar  relationship  exists  between  Eqs.  (123 )  and 
(125),  except  for  sign.  This  property  is  common  to  all  finite  difference 
analyses  of  the  cantilever  beam  and  results  from  the  complementary  character 
of  y  and  M  and  also  of  9  and  V. 

Setting  h  =  1,  one  obtains  the  characteristic  equations  for  the  system: 


*10  -  29-01745A8  +  221.62314*6  .  462.52608*4  +  148.71888*2  -  1.37757  =  0 

(126) 

giving  the  normalized  frequencies 

6>i  =  5.51687 

c as  =  22.1688 

0)3  =  59-5564 

0)4  =  99.8114 

0*5  =  155-687 


Note  that  all  except  the  lowest  mode  frequency  are  roughly  equal  to  cor¬ 
responding  frequencies  for  the  6-cell  built-in  case.  An  even  closer  agree¬ 
ment  occurs  in  the  percentage  error  for  corresponding  cases.  In  fact,  for 
large  N  and  comparing  the  higher  modes,  one  finds  essentially  identical  re¬ 
sults  for  the  built-in  and  cantilever  beams.  (See  Tables  II  and  III.) 
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C.  ORTHOGONALITY 

We  have  seen  that  the  use  of  approximations  of  order  h2  leads  to  orthog¬ 
onal  modes.  Also,  one  can  obtain  a  lumped  physical  system  composed  of  linear, 
bilateral  elements  that  is  represented  exactly  by  the  equations  of  motion. 

This  is  no  longer  true  for  the  h4  approximation  using  one-sided  differences 
at  the  boundaries.  If  one  writes  the  equations  of  motion,  the  m  matrix  is 
symmetric  but  the  k  matrix  is  not.  [See  Eq.  (120),  for  example.]  Thus  the 
system  is  nonphysical,  i.e.,  a  passive  system  of  linear  bilateral  elements 
cannot  be  obtained  for  which  the  equations  are  an  exact  mathematical  repre¬ 
sentation.  Also,  the  modes  are  not  orthogonal,  implying  that  one  cannot  ob¬ 
tain  a  coordinate  transformation  that  diagonalizes  the  m  and  k  matrices 
s imultaneous ly . 

Tb  illustrate  this  point,  consider  again  the  6-cell  built-in  beam  for 
the  case  of  even  symmetry.  For  each  mode,  one  can  calculate  an  amplitude 
ratio  that  indicates  the  relative  magnitude  of  the  motions  of  y2  and  y3. 


Aa  _  5 .22099-A2 
Aa  2.77553 


(127) 


The  results  for  this  case  are  given  by  the  matrix 


1.00000 
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1.00000 

-0.88897 


(128) 


where  each  column  corresponds  to  a  natural  mode  and  the  amplitude  of  the  y2 
motion  is  arbitrarily  set  equal  to  unity. 

The  original  mass  matrix  for  this  system  is 
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The  deviations  from  a  unit  matrix  in  this  case  are  brought  about  by  the  elimina¬ 
tion  of  yj,  using  the  constraint  equation 


yi  =  | ya  ■  k 


Hie  stiffness  matrix  for  the  original  system  is 


Consider  now  the  coordinate  transformation 


(130) 


(131) 


(y)  =  [A]{y] 


(132) 


Hie  mass  matrix  in  terms  of  the  barred  coordinate  system  is 


[m]  *  [A]T[m][A] 


4.05185  -0.50770 

-0.50770  I.8567I 


and  the  stiffness  matrix  is 


[E] 


[AjT[k][A] 


I.58547  -5.90542 

-0.19842  14.27502 


(133) 


(134) 


Hius  we  see  that  neither  the  mass  nor  the  stiffness  matrix  is  diagonalized 
by  the  transformation  to  generalized  modal  coordinates,  and  therefore  the  modes 
are  not  orthogonal .  Furthermore ,  the  k  matrix  is  not  symmetric . 

We  have  defined  orthogonality  of  modes  in  terms  of  a  simultaneous  dlagonaliza- 
tion  of  the  mass  and  stiffness  matrices  by  a  coordinate  transformation.  It  should 
be  noted,  however,  that  a  nonorthogonality  of  modes  in  this  sense  does  not  nec¬ 
essarily  imply  dependence  or  coupling  of  modes.  In  the  above  case,  for  example, 
the  matrix  [BT]”1  [It]  turns  out  to  be  diagonal. 
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VI.  THE  USE  OF  SYMMETRY  IN  REPRESENTING  BOUNDARY  CONDITIONS 

The  use  of  an  h4  approximation  in  a  straightforward  manner  has  been  shown 
to  require  the  use  of  one-sided  differences  near  the  ends,  whereas  central 
differences  can  be  used  in  middle  portion.  Another  approach  to  the  problem 
of  representing  boundary  conditions  is  to  assume  a  virtual  or  image  beam  ex¬ 
tending  beyond  the  actual  beam  and  to  use  values  of  y,  9,  M,  and  V  in  this 
imaginary  beam,  as  necessary,  in  the  standard  central  difference  equations . 
This  approach  will  now  be  investigated. 


A.  SYMMETRY  ASSUMPTIONS 

In  establishing  the  basic  approach  to  the  problem,  we  will  again  use  the 
central  difference  approximation 
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using  staggered  stations  such  that  y  and  M  are  calculated  at  integer  stations 
whereas  9  and  V  are  calculated  at  half-integer  stations .  The  beam  ends  occur 
at  half-integer  stations. 


Let  us  now  consider  appropriate  symmetry  assumptions  for  a  cantilever  beam, 
chosen  to  illustrate  free  and  built-in  end  conditions.  (See  Fig.  10.)  First 
we  note  that  the  assumption  of  even  symmetry  in  y  about  the  built-in  end  will 
give  zero  slope  at  that  end,  using  Eq.  (99)*  (The  assumption  of  odd  symmetry 
here  would  not  give  the  proper  result.)  Successive  differentiation  of  y  with 
respect  to  x  gives  alternating  odd  and  even  functions,  as  shown  by  the  dashed 
lines  where  they  deviate  from  the  solid  lines.  Thus  the  assumed  beam  shows 
a  discontinuity  in  V  (and  the  slope  of  M).  However,  one  can  see  from  the  basic 
approximation  given  in  Eq.  (99)  that  the  presence  of  discontinuities  will 
amount  to  an  increase  in  the  magnitude  of  the  higher  derivatives  at  this  point 
and  thereby  increase  the  error  in  the  approximation. 

To  avoid  these  errors  due  to  discontinuities,  we  will  arbitrarily  assume 
continuous  curves  at  the  boundaries  even  though  these  assumptions  are  not 
physically  consistent.  We  are  interested  in  obtaining  the  best  possible 
estimates  for  the  derivatives  of  y  within  the  actual  beam,  and  therefore  we 
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make  the  most  promising  assumptions  concerning  the  virtual  beam. 

The  symmetry  assumptions  that  we  will  use  can  be  summarized  as  follows: 

(1)  If  the  derivative  of  the  function  approaches  zero  as  the  end  of  the 
actual  beam  is  approached,  assume  even  symmetry  about  the  end. 

(2)  If  the  derivative  of  the  function  does  not  approach  zero  as  the  end 
is  approached,  assume  odd  symmetry  about  the  end,  possibly  with  an  offset  to 
avoid  discontinuities. 

Ihe  procedure  by  which  these  assumptions  are  incorporated  into  the  equa¬ 
tions  of  motion  is  quite  straightforward  except  for  the  extrapolation  of  M 
at  a  built-in  end  or  the  similar  extrapolation  of  y  at  a  free  end.  Taking 
the  case  of  the  built-in  end,  a  functional  form  must  be  assumed  for  M  near 
the  boundary.  Referring  to  the  analytical  solution  near  a  built-in  end,  as 
given  by  Eq.  (22)  or  (49),  we  find  that  the  magnitude  of  M  varies  according 
to  a  linear  plus  a  fourth-order  term  in  x  plus  other  higher-order  terms .  So 
a  logical  choice  of  the  functional  form  of  M  near  the  boundary  might  be 


M  =  a0  +  ExX  +  e,4x4 


(135) 


where  the  last  term  takes  the  plus  or  minus  sign  depending  on  whether  x  is 
positive  or  negative.  This  assumption  results  in  the  equations 


Mo  =  ~  [503Mi  -  310Ma  +  39M3]  (156) 


and 

M-1  =  252  [T55Ml  "  542Ms  +  59Ma]  U37) 

where  the  built-in  end  occurs  at  station  £. 

Another  possibility  is  to  assume  that  M  is  an  offset  odd  function  having 
a  linear  plus  cubic  variation  with  x.  Ihus  we  might  assume  that 


M 


b0  +  bix  +  b3x3 


(138) 
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j.uib  assumption  avoids  having  discontinuous  derivatives  at  x  ■  0  and  gives 
numerically  simpler  equations.  In  this  case 


Mo  =  l[7Mi.-  5 Mb  +  M3] 


(139) 


M.i  =  l[10Mi  -  8Ma  +  M3] 


(140) 


Computations  were  run  for  many  cases,  using  either  Eq.  (135)  or  (138),  and 
the  results  were  compared.  The  differences  were  fairly  small  and,  if  anything, 
favored  the  linear  plus  cubic  assumption.  All  results  given  in  this  report, 
assuming  symmetry  at  the  boundaries,  are  based  upon  a  linear  plus  cubic  varia¬ 
tion  of  M  near  a  built-in  end,  i.e.,  upon  Eqs.  (139)  and  (l4o). 

3he  assumptions  concerning  the  form  of  y  near  a  built-in  end  are  the 
same  as  we  used  previously  for  the  h4  case.  Again  we  can  eliminate  yi  by 
using  the  equation 


yi  =  § ys  ■  k ys 


(130) 


Of  course  a  similar  relation  exists  for  the  value  of  M  at  the  station 
nearest  a  free  end. 


B.  EQUATIONS  FOR  A  BUILT-IN  BEAM 

1.  General  Case.— We  have  seen  how  we  can  use  symmetry  to  extrapolate 
the  variables  of  the  problem  beyond  the  actual  limits  of  the  beam.  Using  these 
extrapolated  values  and  the  basic  centred,  difference  approximation  of  Eq.  (99) > 
we  caua  write  genered.  equations  for  a  built-in  beam  simileur  to  those  given  by 
h4  approximation  in  Eqs .  ( 106 )  to  ( 109 ) •  Ihey  are 
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A  comparison  of  the  above  equations  with  those  for  the  h4  approximation 
as  given  by  Eqs.  (106)  to  (109)  shows  considerable  similarity.  It  should  be 
noted,  however,  that  the  use  of  symmetry  at  the  boundary  never  requires  more 
than  four  inputs  per  amplifier  in  the  analog  computer  circuit,  whereas  the  h4 
approximation  requires  up  to  five  inputs  per  amplifier. 

2.  Ihe  6-Cell  Case.— As  an  example  of  a  mode  calculation  for  a  specific 
case,  consider  again  the  modes  with  even  symmetry  about  the  center  for  a  6- 
cell  beam  built-in  at  both  ends.  The  equations  of  motion  (using  symmetry 
assumptions  at  the  boundaries )  are 


■fr  CNVO  O 


The  roots  are 


A1;2  =  0.62828,  3.25078 


and  the  normalized  roots  are 


0>!  =  22.6182 

Ofe  =  117.028 


Recalling  that  the  corresponding  normalized  frequencies  for  the  con¬ 
tinuous  beam  are 


Pi  =  22.3733 

fe  =  120.903 


we  see  that  the  above  results  give  a  somewhat  larger  error  for  the  first  mode 
than  does  the  h4  approximation,  but  the  error  in  the  higher  mode  is  much  smaller. 

An  alternate  method  can  be  used,  particularly  with  the  use  of  a  digital 
computer,  to  set  up  the  equations  of  motion  and  to  solve  them.  Ihis  method, 
which  indicates  more  clearly  the  symmetry  assumptions  in  its  formulation,  is 
given  in  Appendix  A. 

3.  Results . — Hie  results  of  calculations  for  built-in  and  cantilever  beams, 
using  symmetry  assumptions  at  the  boundaries,  are  summarized  in  Tables  II  and 
III  and  in  Figs.  17  and  18. 

Comparing  these  results  with  those  for  the  h4  approximation,  we  note  that 
for  a  practical  range  of  N  (say  6  to  16)  the  use  of  symmetry  gives  lower  ac¬ 
curacy  on  the  low  frequency  modes  but  better  accuracy  on  the  higher  modes.  For 
example,  consider  the  case  of  an  8-cell  cantilever  beam  having  a  total  of 
seven  modes.  Using  symmetry  at  the  boundaries,  the  first  four  modes  have  a 
frequency  error  well  within  1$  and  even  the  sixth  mode  has  an  error  of  only  10$, 
approximately.  By  contrast,  the  use  of  h4  approximations  at  the  boundaries 
gives  better  accuracy  for  the  first  three  modes  but  the  fourth  mode  has  an 
error  of  nearly  5$  and  the  sixth  mode  has  an  error  of  approximately  20$. 


c.  ORTHOGONALITY 


It  turns  out  when  boundary  conditions  are  imposed  using  symmetry  that 
the  resulting  modes  of  vibration  are  nearly,  but  not  exactly,  orthogonal. 
Ihis  means  that  this  approach  also  results  in  a  nonphysical  system  in  the 
sense  that  a  lumped  model  cannot  be  built  of  linear  bilateral  elements. 


To  compare  the  orthogonality  properties  with  those  for  the  h4  approxi¬ 
mation,  let  us  consider  again  the  example  of  the  6-cell,  built-in  beam. 

The  mass  matrix  is  again 


[m] 


Hie  stiffness  matrix  is 
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-2 
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225 
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(129) 


[k]  «  [m][K] 


8.48426  -4.63270 

-4.47131  2.95823 


(15D 


Hie  amplitude  ratios  are  calculated  from 


Aa  _  8.Q4779-A2 
Ax  4.39000 


(152) 


yielding  the  modal  matrix 


[A] 


1.00000  1.00000 
1.7*329  -0.57399 


(153) 


In  the  modal  coordinate  system  the  mass  matrix  is 
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[nr]  -  [A]T[m][A] 


4.06232  0.Q3676 

0.03676  1.38957 


The  stiffness  matrix  is 

PH  =  [AjT[k  ]  [A] 


I.60356  O.38348 

0.01451  14.68446 


(154) 


(155) 


Thus  neither  matrix  is  diagonalized,  by  the  transformation,  but  the  result 
is  considerably  closer  than  was  the  case  for  the  h4  approximation,  as  can  be 
seen  from  Eqs.  (135)  and  (134). 
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VII.  THE  PASSIVE  CIRCUIT  METHOD 


We  have  seen  in  the  case  of  the  simple  approximation  of  order  h2  that  an 
electrical  circuit  exists  which  the  difference  equations  describe  exactly. 
However,  when  one  goes  to  higher-order  differences  and  :ses  the  methods  we 
have  described  for  establishing  the  boundary  conditions,  one  finds  that  the 
resulting  equations  are  nonphysical,  i.e.,  no  passive  linear  circuit  analogy 
exists. 

It  can  be  shown  that  the  nonphysical  character  of  the  higher-order  methods 
discussed  thus  far  is  due  to  the  boundary-condition  representation  rather  than 
the  basic  difference  approximation.  In  this  section  a  passive  circuit  analogy 
will  be  obtained  for  the  basic  higher-order  approximations  of  Eqs.  (100)  to 
(103).  By  imposing  boundary  conditions  on  this  circuit  and  writing  the  cor¬ 
responding  equations,  one  can  obtain  a  set  of  mutually  orthogonal  modes. 


A.  PASSIVE  CIRCUIT  ANALOGY 

We  have  been  using  the  following  finite  difference  approximations  to  the 
beam  equations. 


9n-*i 

_  1 
24h 

y^n-i  “  27yn  +  2Tyn+1  -  yn+2J 

(100) 

M  = 
n 

2kh  \ 

'n-j  '  2Ten-J  +  279n+J  '  Sn»|) 

(101) 

Vn+£ 

=  1 
24h 

-  27Mn  +  27M,»i  -  *W) 

(102) 

yn  = 

—  (■ 

24h  \ 

'Vj  *  27Vj  •  27V*  *  O 

(103) 

Setting  h  =  1,  and  representing  generalized  velocities  by  voltages  and  gen¬ 
eralized  forces  by  currents,  one  can  show  that  the  passive  circuit  of  Fig.  11 
is  described  by  these  equations.  Note  that  three  transformers,  one  inductor, 
and  one  capacitor  sure  required  per  cell. 

The  rather  large  transformer  requirement  of  Fig.  11  compared  with  the  h2 
approximation  of  Fig.  3(t>)  makes  it  impracticsLL  to  use  this  passive  circuit 
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directly.  Instead  we  can  write  the  circuit  equations  and  solve  them  vising  a 
differential  analyzer  or  a  digital  computer.  The  equipment  requirements  in 
this  case  are  almost  the  some  as  for  the  other  higher-order  methods  but  are 
actually  slightly  less  due  to  the  details  of  boundary  representation. 


B.  THE  BUILT-IN  BEAM 

1.  Symmetry  Assumptions . — To  obtain  a  proper  representation  of  a  built- 
in  end,  it  is  convenient  to  make  certain  symmetry  assumptions  consistent 
with  the  given  boundary  conditions.  We  will  assume  that  the  deflection  curve 
has  even  symmetry  about  the  end  and  zero  slope  at  the  boundary.  Thus  y  and 
M  are  even  functions  of  x,  while  6  and  V  are  odd  functions  of  x,  as  shown 
by  the  dashed  lines  of  Fig.  10. 

The  circuit  that  has  been  used  to  represent  a  built-in  end  is  shown  in 
Fig.  12.  The  beam  end  is  at  the  center  of  the  diagram  (station  £)  with  the 
virtual  beam  on  the  left  and  the  actual  beam  on  the  right.  The  autotrans¬ 
former  at  the  bottom  of  the  figure  represents  the  constraint  relating  y1 
and  y2- 

2.  Constraint  Equation. — It  can  be  shown  for  a  uniform  beam  that 

V*  '  *  9,“  *  '  r"+J)  *  fS  S 

So  if  we  assume  that  y  has  even  symmetry  about  the  built-in  end  at  station  1 
and  also  set  y^  equal  to  zero,  we  obtain  the  approximate  equality 

yi  =  |  ya  (157) 

This  is  used  as  the  constraint  equation  in  the  passive  circuit  approach. 

It  can  be  seen  from  Fig.  12  that  the  autotransformer  representing  the  con¬ 
straint  equation  causes  a  certain  relationship  to  exist  between  voltages  and 
also  between  currents.  In  mechanical  terms,  a  constraint  equation  relating 
coordinates  also  implies  forces  of  constraint  to  be  exerted  on  the  system  in 
such  a  way  that  no  work  is  done.  These  forces  of  constraint  have  been  ignored 
in  the  previous  two  approaches  to  the  problem  and  this  resulted  in  a  lack  of 
symmetry  in  the  stiffness  matrices. 


n+2 


(156) 
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To  calculate  the  effect  of  these  forces  of  constraint,  consider  a  beam 
described  by  the  matrix  equation 


[m]{y}  -  p](V)  (158) 

Nov  suppose  a  constraint  is  imposed  as  described  by 

yi  -  (c)(y)  (159) 


where  the  summation  is  over  all  y  except  yi.  (Hhe  parentheses  denote  a  row 
matrix. )  Let  us  include  the  constraint  forces  in  the  analysis  and  also  separate 
the  yx  equation  from  the  others.  Ihen  we  obtain 


mm  =  (Di)lV)  +  Fx 

(160) 

MW)  =  [D](V)  +  (F) 

(161) 

where  the  F  are  constraint  forces  acting  on  the  beam.  Note  that  Eq.  (l6l)  does 
not  contain  yx  and  the  D  matrix  does  not  contain  a  Dx  row. 

T 

Now  multiply  (160)  by  yx  and  premultiply  (l6l)  by  (y)  .  Adding,  we  obtain 
yxmiyi  +  (y}T[m]{y)  =  yi(Di)(V)  +  {y}T[D](v}  +  Fiyi  +  (y}T(F)  (162) 


Since  the  constraint  does  no  work,  we  can  write 


Fiyi  +  (y}T(F)  =  0 


(163) 


and  the  last  two  terms  of  (162)  drop  out.  Finally,  substituting  for  yi  from 
Eq.  (159),  dividing  out  (y } T,  said  rearranging,  we  obtain 
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n  owo  -4- 


mi(C)(y}(C)T  +  Cm](y)  -  (D!){V}(C)T  +  [D](VJ 
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(164) 
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This  gives  the  equations  of  motion  with  yi  omitted  and  constraint  forces  in¬ 
cluded.  For  example,  in  the  present  case  with  a  simple  constraint  given  by 
Eq.  (157),  the  effective  mass  at  yz  is  82/81  times  its  value  before  the  con¬ 
straint  was  applied.  More  complicated  constraints  will  result  in  inertial 
coupling,  but  the  effective  inertia  matrix  will  always  be  symmetric. 

3.  General  Case. — The  complete  equations  of  motion  for  a  built-in  beam 
can  be  written  now  using  the  passive  circuit  method.  They  are 
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The  amount  of  analog  computer  equipment  needed  for  the  passive  circuit 
method  is  slightly  less  than  for  either  of  the  other  higher-order  methods. 
There  are  two  reasons  for  this.  First,  the  value  of  V  at  a  "built-in  end  (or 
of  0  at  a  free  end)  need  not  be  calculated.  Second,  the  signs  are  optimum  in 
that  no  extra  amplifiers  are  required  for  inversion  at  any  point. 

The  equations  of  motion  for  stations  near  a  free  end  can  be  obtained  most 
easily  from  the  above  equations  by  exchanging  the  roles  of  y  and  M  and  also 
6  and  V.  Of  course  they  could  also  be  obtained  directly  from  the  passive  cir¬ 
cuit,  remembering  the  constraint  equation  relating  the  values  of  M  at  the  last 
two  stations . 


4.  The  6 -Cell  Case. — Again  let  us  consider  the  even  symmetry  modes  of 
a  6-cell  beam  built-in  at  both  ends.  Using  the  passive  circuit  approach,  the 
equations  of  motion  are 
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Setting  h  =  1,  we  can  solve  for  the  eigenvalues  from 


(8.4T527-X2)  -4.54772 

-4.60386  (2.99140-X2) 


0 


or 


(172) 


(175) 


X4  -  11.46667X2  +  4.41583  =  0 


(174) 


The  normalized,  roots  are 


0>!  =  22.7395 

cite  3  119-765 

In  this  case  the  first  mode  frequency  is  1.6$  high  and  the  second  mode  fre¬ 
quency  is  0.94$  low. 


C.  RESULTS 

Die  results  of  mode  frequency  calculations  for  a  uniform  beam  using  the 
passive  circuit  method  are  summarized  in  Tables  II  and  III  and  in  Figs.  19  and 
20. 


Comparing  these  results  with  those  obtained  previously,  using  other  methods, 
we  note  that  the  passive  circuit  method  gives  very  good  over-all  accuracy  for 
cases  where  the  number  of  cells  is  9  or  more.  On  the  other  hand,  the  symmetry 
method  appears  to  give  generally  better  results  for  smaller  values  of  N.  Of 
course,  the  h4  method  is  still  best  for  accuracies  of  the  order  of  0.25$  or 
better. 
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One  way  of  comparing  the  various  methods  Is  to  obtain  an  approximate  number 
of  cells  required  per  wavelength  to  obtain  a  given  accuracy.  For  1$  error, 
the  ordinary  h2  approximation  requires  12  cells  per  wavelength.  In  contrast, 
the  h4  approximation  requires  6-l/2  cells  per  wavelength,  the  symmetry  meth¬ 
od  requires  5  cells  per  wavelength,  and  the  passive  circuit  method  requires 
4-1/2  cells  per  wavelength. 


D.  ORTHOGONALITY 

In  contrast  with  the  higher-order  methods  considered  previously,  the  pas¬ 
sive  circuit  method  produces  orthogonal  modes.  The  orthogonality  property  is 
a  result  of  the  symmetry  of  the  original  m  and  k  matrices. 

We  can  illustrate  this  point  by  referring  again  to  the  6-cell  case.  From 
Eq.  (173)  we  see  that  the  amplitude  ratios  for  the  two  modes  can  be  calculated 
using 


Afe  8.47527-A2 
Ai  =  4.54772 


(175) 


The  modal  matrix  is 


* 


[A] 


1.00000  1.00000 

1.77590  -0.57005 


(176) 


The  original  mass  matrix  is 


The  corresponding  stiffness  matrix  is 


[k] 


8-57990  -4.60386 

-4.60386 


(177) 


2.99140 


(178) 
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After  the  transformation 

(y)  -  [A)(y}  (132) 


Hie  diagonal  nature  of  these  matrices  implies  orthogonality  of  modes,  and  also 
that  the  motion  of  each  of  the  modal  coordinates  is  independent  of  the  others. 

It  should  be  noted  in  passing  that  the  matrix 


[A2  ]  =  [AfVHA] 


(181) 


is  diagonal  for  all  three  methods  and  the  numbers  along  the  main  diagonal  sire 
the  squares  of  the  mode  frequencies.^-®  Ihis  does  not  imply  orthogonality, 
however,  as  we  have  defined  it. 
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VIII .  SOURCES  OF  ERROR 


In  this  section  we  will  consider  the  principal  sources  of  error  in  the 
methods  that  have  been  suggested  for  the  solution  of  beam  vibration  problems. 
These  errors  are  (1)  finite  difference  errors  due  to  the  basic  difference  approx' 
imation;  (2)  finite  difference  errors  in  the  representation  of  boundary  con¬ 
ditions,  and  (5)  component  errors.  The  error  curves  that  have  been  presented 
so  far  represent  a  combination  of  error  sources  (l)  and  (2)  above  for  various 
specific  cases.  Now  we  will  consider  them  with  somewhat  more  generality. 


A.  FINITE  DIFFERENCE  ERRORS 

1.  The  h2  Approximation .  We  have  seen  previously  for  a  pinned  beam  that, 
when  proper  symmetry  assumptions  are  made,  there  is  no  additional  error  due  to 
the  method  of  representing  boundary  conditions  and  so  the  entire  computed 
error  is  due  to  the  basic  finite  difference  approximation.  So  consider  now 
a  uniform  beam  pinned  at  each  end  at  an  integer  station. 

Recall  again  the  basic  h2  approximation 


(60) 
n 

Repeated  application  of  Eq.  (60)  yields 


dy  ^  ~yn  +  yn+i  h2  d3y 

dx  i  h  24  dx3 

n+? 


df£ 

dx2 


yn-i  '  2yn  *  yn+i  hf  d^ 
h2  12  dx4 


(65) 


Similarly,  one  can  obtain 


yn-g  -  4yn-i  +  6yn  -  4yn+!  +  yn+g  hf.  d^ 
h4  "6  ax6 


(182) 


which  is  the  result  of  applying  Eq.  (60)  four  times  in  sequence 


on  vo  t» 
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We  will  again  consider  the  equation  for  the  continuous  beam  in  simplified 

form. 


0  (6) 

ax*  at2 


As  we  have  seen  previously,  the  solution  for  the  continuous  pinned  beam  is  of 
the  form 


y  =  Z  sin n /ftT  x  (G.  cos  f^t  +  H.  sin  0^t)  (183) 

i 

where 


^  =  i2n2  (i  =  1,  2,  5,  ...)  (55) 

For  the  case  of  the  ith  mode,  we  can  substitute  Eqs.  (182)  and  (185)  into  (6), 
obtaining 

7  (yn-2  -  K-i  +  6yn  -  4yn+i  +  w)  +  f  P?yn  -  ^yn  =  0 

or 

7  (yn-2  '  4yn-i  +  6yn  "  4yn+1  +  yn+2)  =  ‘  f  * )  <l84) 

Now  the  left-hand  side  of  Eq.  (184 )  is  of  the  same  form  as  the  approximation 
to  d4y  used  in  the  finite  difference  solution,  therefore,  the  corresponding 
dx4  n 

finite  difference  equation  is 

~  (yn-2  -  4yn-i +  6yn  -  4yn+i +  yn«)  =  4?*. 


(185) 
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where  is  the  finite  difference  eigenvalue  corresponding  to  S^.  So  the  right 
sides  of  Eqs.  (l84)  and  (185 )  can  be  equated  if  sinusoidal  mode  shapes  are 
assumed  in  each  case.  Finally,  noting  that 


h  *  i 


(186) 


we  obtain 


'  (1">4  £  -^r)* + 


or 


=  (iJt)“ 


[-my*-] 


(187) 


where  i  is  the  mode  number. 

2.  The  h4  Approximation . — A  similar  approach  can  be  used  to  compute  the 
frequency  error  for  a  pinned  beam  using  higher-order  differences.  The  basic 
approximation  is 


k 

dx 


n+A 


(99) 


n+J 


Using  successive  repetitions  of  Eq.  (99) >  we  obtain 


£ 1 

8x2 


n  =  (Ikhj*  (yn-3‘  5^n-2+785yn-x-li+6oyn+  g  g 


(188) 


and 


*2 

8x4 


(24h) 


+  . . .  +  y  *^  + 
e  •'n +6/ 


160  8x8 


(189) 
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Proceeding  in  a  manner  similar  to  that  for  the  h2  approximation,  ve  ob¬ 


tain 


or 


0)7 


A 


. 

160 


A 


(190) 


Thus  we  see  that  the  basic  finite  difference  approximation  tends  to  give 
frequencies  that  are  too  low.  The  first  error  term  in  Eq.  (190)  gives  an 
accurate  estimate  of  the  error  only  for  small  i/N  ratios  tod  more  terms  are 
needed  for  most  practical  cases.  The  actual  error  curves  are  plotted  in  Fig. 
21  for  the  h4  approximation. 

3.  Maximum  Errors . — From  the  tabulated  results  of  the  calculations  of 
mode  frequencies,  it  is  seen  that,  as  the  i/N  ratio  increases,  the  error  also 
increases.  In  the  limit  as  i/N  approaches  unity,  the  error  is  independent 
of  the  boundary  conditions  and  depends  only  on  the  basic  approximation.  This 
limiting  value  of  the  error  could  be  obtained  by  evaluating  the  complete 
series  indicated  in  Eq.  (l87)  or  (190)  for  i/N  «  1.  A  simpler  way  to  obtain 
this  result  is  to  note  that  as  i/N  approaches  unity  the  motion  at  adjacent 
stations  is  equal  in  amplitude  but  opposite  in  sign.  Knowing  this,  one  can 
calculate  the  frequency  directly. 

For  example,  using  the  h2  approximation  and  assuming  h  >1,  we  find  that 


yn  +  i6yn  -  0 


(191) 


yielding  the  frequency 
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Hie  corresponding  frequency  for  the  continuous  beam  is  the  same  as  for  a  pinned 
beam  of  unit  length  vibrating  in  its  fundamental  mode’,  namely,  nz .  So  the 
maximum  error  fcr  the  h2  approximation  is 

W  -  m  S’-1*72*  <192) 

Jt 

Performing  a  similar  calculation  for  the  h4  case,  we  obtain 

yn  *  (jj  ■  0  (1») 

or 


A 


hi 

9 


giving  a  maximum  error 


€max 


9wg-49 

9n2 


44.836$ 


(194) 


B.  BOUNDARY-CONDITION  ERRORS 

We  have  seen  that  an  analysis  of  uniform  pinned  beams  with  various  numbers 
of  cells  gives  directly  the  finite  difference  error  due  to  the  basic  centred 
difference  approximation.  Furthermore,  we  note  that  the  eigenvalues  of  the 
pinned  case  alternate  with  the  eigenvalues  of  either  the  cantilever  or  built- 
in  case.  Consequently,  if  the  method  of  applying  boundary  conditions  intro¬ 
duced  no  additional  error,  one  would  expect  the  error  curves  for  the  cantilever 
or  built-in  beam  to  run  between  and  roughly  parallel  to  the  error  curves  for 
the  pinned  beam.  Instead,  one  finds  that  all  the  boundary-condition  methods 
considered  here  result  in  raising  the  natural  frequencies.  His  causes  a  gen¬ 
eral  improvement  in  accuracy  for  i/N  values  greater  than  about  0-3  for  the 
symmetry  and  passive  circuit  approximations.  For  the  h4  method  the  accuracy 
is  improved  for  all  i/N. 
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For  the  lower  modes  and  large  N,  the  error  for  the  three  built-in  or 
cantilever  ewes  mentioned  above  exceeds  the  error  for  the  corresponding  pinned 
case  because  the  boundary-condition  errors  are  of  lower  order  in  h  than  the 
basic  central  difference  errors,  and  ultimately  must  predominate.  For  example, 
in  the  equations  for  a  built-in  beam  using  the  simple  "h2 "  approximation,  we 
find  the  expression 


6a  =  —  (195) 

2  h 

[See  Eq.  (82).]  The  error  in  this  case  is  actually  of  order  h,  as  can  be  ascer¬ 
tained  by  the  Taylor  expansion  method. 

Similarly,  the  passive  circuit  and  symmetry  methods  of  approximating  a  built- 
in  end  use  the  expression 


which  has  a  first  error  term  of  order  h.  On  the  other  hand,  each  equation  used 
in  the  h4  approximation  is  truly  of  order  h4,  and  therefore  the  boundary- 
condition  errors  do  not  tend  to  dominate  for  large  N. 

A  study  of  the  error  curves  presented  in  Figs.  l6  and  21  It* ads  one  to  the 
conclusion  that  the  first  error  term  alone  does  not  give  a  reliable  estimate 
of  the  actual  errors  encountered.  The  i/N  ratio  must  be  the  order  of  0.1  or 
less  before  the  first  error  term  clearly  predominates  and  by  this  time  roundoff 
or  component  errors  mask  the  desired  truncation  error. 

To  illustrate  this  point,  it  is  of  interest  to  calculate  the  rate  at 
which  the  actual  truncation  errors  decrease  with  increasing  N  as  N  approaches 
16,  i.e.,  at  the  right-hand  edge  of  the  error  diagram.  Rough  calculations 
indicate  that  the  h2  and  symmetry  methods  have  an  error  that  is  decreasing 
as  h1'8,  while  the  error  for  the  passive  circuit  method  is  decreasing  as  h3'0 
and  that  for  the  h4  method  is  decreasing  as  h4,7  for  increasing  N. 

The  conclusion  from  the  preceding  analysis  is  that  the  method  of  repre¬ 
senting  boundary  conditions  should  be  chosen  on  the  basis  of  the  require¬ 
ments  of  the  problem  and  actual  error  curves ,  rather  than  assuming  that  an 
h4  approximation  is  necessarily  more  accurate  than  an  h3  or  lower  approxima¬ 
tion  over  the  practical  range  of  N.  Our  results  show  that  so-called  higher- 
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order  methods  are  generally  superior  to  simpler  approximation  methods.  How¬ 
ever,  the  choice  between  the  various  methods  of  representing  boundary  condi¬ 
tions  should  be  made  on  the  basis  of  the  problem  at  hand. 


C.  COMPONENT  ERRORS 

In  this  section  we  will  consider  the  sensitivity  of  the  calculated  roots 
or  frequencies  to  small  shifts  or  inaccuracies  of  the  component  values.  Ibis 
will  be  accomplished  by  perturbing  one  of  the  matrix  elements  and  calculating 
the  resulting  frequency  shift. 

Suppose  we  write  the  characteristic  equation  in  the  form 


S(tj,  bij)  =  0 


(197) 


where  the  eigenvalue,  tj  =  A2,  is  the  square  of  the  frequency  (rad/sec)  of  one 
of  the  modes  and  b^j  is  the  matrix  element  that  is  to  be  perturbed.  For 
small  deviations  about  the  reference  values  bj_  j0  we  can  expand  in  the 
series 

Ab^  +  ...  (198) 

o 

where 

tj  *  +  Aq 

*ij  =  ^i  Jo  + 


S(q,  bij)  -  S(t|0,  bijQ)  + 


dS 


Atj  + 


as 

^biJ 


From  Eqs.  (197)  and  (198)  we  obtain  an  expression  for  the  shift  in  the  root. 


Atj 


bS 

1 

dbit 

L 

as 

dn 

0 

U 


(199) 


To  illustrate  this  method,  let  us  make  an  error-sensitivity  computation 
for  a  6-cell  beam  that  is  built-in  at  both  ends  using  symmetry  assumptions  at 
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the  ends  and  assuming  even  symmetry  about  the  center.  We  have  seen  that  the 
equations  of  motion  can  he  written  In  the  form 


m  +  [K](y)  -  o 


(79) 


where 


[K]  =  -  [B^tBaHBaHB!] 


(80) 


For  the  case  we  are  considering,  the  characteristic  equation  is  of  the  form 


(Kn-n)  K!2 

Kai  (Kj.i-tj) 


0 


or 


S  =  T|2  -  (Kh  +  Kgg)t|  +  (KnlQg2  -  KtpK pi )  *  0 


(200) 


Let  us  suppose  that  an  element  in  the  Bi  matrix  is  to  he  varied.  We  can  write 


[K]  =  -  [T][BX]  (201) 

where 

[T]  -  [B4][B3][B2]  (202) 


Consider  now  the  case  at  hand  where 
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4775  9 


-6025  5841J 


and  evaluate  the  sensitivity  to  a  perturbation  of  b21.  Evaluation  of  Atj  ac¬ 
cording  to  Eq.  (199) >  assuming  h  =  1,  gives 

Atji  -  0.157010  Ab21 
At)2  =  4.22168  Ab2i 

for  the  first  and  second  modes,  respectively.  In  this  case  the  actual  roots 
are 


(Bil 


5400h 


ru  =  0.394741 

q2  =  IO.5676 


So  the  fractional  changes  in  the  roots  are 

£ll  =  _  0.443791  ^ZX 

^21 

=  .  0.445731  £3221 

b21 

This  implies  that  a  1$  increase  in  the  magnitude  of  b2i  will  cause  a  0.222$ 
decrease  in  the  first  mode  frequency  and  a  0.223$  decrease  in  the  second  mode 
frequency. 

8  Q 

It  has  been  pointed  out  by  Clymer,  Fisher,7  and  others  that,  in  general, 
the  repeated  use  of  finite  difference  approximations  to  the  first  spatial 
derivative  will  result  in  less  component  error  sensitivity  than  using  approx¬ 
imations  to  the  second  or  fourth  derivatives.  To  check  this  point,  let  us 
calculate  the  error  sensitivity  for  the  same  problem  using  approximations  to 
the  second  derivative.  Splitting  the  matrix  multiplication  differently,  we 
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[K]  = 

-  [T][B] 

(203) 

where 

[T]  = 

[B4][B3] 

(204) 

[B]  = 

[B2][B!] 

(205) 

In  this  case,  for  h  =  1,  we  find  that 


[B] 


- 

58.1271 

-130.2083 

72.0812 


-2.5037 

70.3125 

-67.8088 


x  10 


-5 


and 


[T] 


-1981  4115  -2134 
152  -2182  2030 


Again  we  will  calculate  the  error  sensitivity  for  t>2i.  Performing  the  numerical 
computations  as  indicated  in  the  equation,  we  obtain 


Atu  *  -77.6587  Ab21 

At)2  =  -4037-3**  Ab21 

corresponding  to 

^  =  O.256163 

Tli 


i 
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-  0.497460  ^81 

^2  bgj 


So  in  this  case  a  1#  increase  in  the  magnitude  of  b2i  will  cause  a  0.128$  in¬ 
crease  in  the  first  mode  frequency  and  a  0.249$  increase  in  the  second  mode 
frequency. 

Comparing  these  results  with  those  for  the  approximation  to  the  first 
derivative,  we  do  not  have  a  very  significant  difference  in  sensitivity.  In 
fact,  neither  case  would  he  particularly  troublesome  from  the  standpoint  of 
component  error  sensitivity.  Further  calculations  of  error  sensitivities  for 
other  matrix  elements  show  comparable  results. 

This  same  case,  using  approximations  to  the  first  derivative,  was  checked 
for  error  sensitivity  on  an  actual  analog  computer  with  results  in  accordance 
with  the  above  calculations.  Percentage  frequency  shifts  were  generally  less 
than  half  the  percentage  component  change. 

As  the  number  of  cells  N  is  increased,  the  differences  in  error  sensitivity 
between  the  first  derivative  and  second  derivative  approximations  should  be¬ 
come  more  apparent.  In  any  case  the  sensitivity  to  component  errors  will  in¬ 
crease  with  N  and  large  sensitivities  should  be  noticed  first  in  the  second 
derivative  approximation. 

As  Fisher  has  pointed  out,  certain  elements  in  the  matrix  will  show  con¬ 
siderably  more  error  sensitivity  than  others .  In  general  the  large  terms  are 
more  sensitive  to  error.  For  our  basic  higher-order  approximation 


dx 


n+* 


—fy  - 

24hV^  n_1 


27y  +  27y 


n 


n+i 


(99) 


we  would  expect  a  larger  sensitivity  to  errors  in  the  coefficients  ±27  than 
in  the  coefficients  ±1  or  even  the  coefficient  24.  The  coefficients  ±1  act  as 
correction  or  trimming  terms  and  do  not  need  to  be  held  to  the  same  percentage 
accuracy. 

The  above  discussion  of  error  sensitivity  applies  to  digital  computers 
when  considering  roundoff  errors.  The  magnitude  of  the  roundoff  error  is,  of 
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course,  dependent  upon  the  number  of  digits  used  in  the  computations  and  the 
cell  number  N.  It  is  also  dependent  upon  the  differencing  scheme  used  and 
the  details  of  the  digital  program.  For  the  methods  used  with  IBM  704  com¬ 
puter  in  calculating  frequencies  for  the  various  cases  given  in  this  report, 
roundoff  errors  of  the  order  of  0 . 01$  were  first  noticed  with  i/N  ratios  of 
approximately  0.1  or  less  on  cantilever  beam  computations.  Because  the  use 
of  symmetry  effectively  halved  the  number  of  degrees  of  freedom,  the  calc  illa¬ 
tions  for  the  built-in  beam  did  not  suffer  appreciably  from  roundoff  errors. 

For  reasons  of  roundoff  error  and  the  general  complexity  of  computation, 

N  ■  16  seems  to  be  a  rough  upper  limit  on  the  number  of  cells  to  be  used  in 
analog  or  digital  computations.  Diis  detracts  somewhat  from  the  utility  of 
the  h4  method  of  representing  boundary  conditions  because  it  is  at  large  N 
that  this  method  is  particularly  advantageous  from  the  standpoint  of  finite 
difference  or  truncation  errors.  On  the  other  hand,  it  accentuates  the  ad¬ 
vantage  of  higher-order  difference  methods  in  general  as  compared  to  the  simpler 
h2  methods,  particularly  in  the  computation  of  the  higher  frequency  modes. 
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IX.  CONCLUSIONS 


The  principal  conclusions  gained  from  the  investigations  reported  here 
are  as  follows. 

1.  The  use  of  higher-order  differences  in  the  solution  of  beam- vibration 
problems  on  the  analog  computer  improves  accuracy  with  essentially  the 
same  amount  of  equipment  as  would  be  required  for  simpler  difference 
methods.  In  the  digital  formulation  of  the  corresponding  eigenvalue 
problem,  a  similar  improvement  in  accuracy  is  obtained  with  essentially 
the  same  computer  effort,  as  measured  by  matrix  sizes  and  numbers  of 
operations.  The  improved  accuracy  is  obtained  by  using  more  inputs 

per  amplifier  in  the  analog  circuit  and  more  nonzero  elements  per  matrix 
in  the  digital  computations. 

2.  The  principal  difficulty  in  the  formulation  of  beam- vibration  problems 
using  higher-order  differences  lies  in  the  proper  representation  of 
boundary  conditions.  This  difficulty  arises  from  the  need  for  obtain¬ 
ing  input  data  from  more  stations  for  each  differencing  operation  to 
obtain  greater  accuracy  than  in  the  simpler  approximations.  As  the 
boundary  is  approached,  some  cf  these  stations  needed  for  the  basic 
higher-order  central  difference  approximation  lie  beyond  the  end  of  the 
beam,  and  so  other  approximations  or  assumptions  must  be  used.  Three 
methods  are  presented  for  the  representation  of  built-in  or  free  ends, 
namely,  (l)  the  use  of  one-sided  differences  of  order  h4,  (2)  the  use 
of  symmetry  assumptions  at  the  boundary,  and  (3)  the  passive  circuit 
method . 

Each  method  has  its  areas  of  special  utility,  as  seen  from  the  re¬ 
sults  of  calculations  for  uniform  beams.  The  h4  method  is  best  suited 
to  relatively  small  i/N  ratios,  corresponding  to  accuracies  of  0.25 $  or 
better.  The  symmetry  method  is  most  advantageous  in  obtaining  moderate¬ 
ly  good  accuracies  (the  order  of  1$  for  all  but  the  highest  modes  with 
a  minimum  of  equipment,  particularly  for  a  cell  number  N  <  9.  The  pas¬ 
sive  circuit  method  has  quite  good  accuracy  in  general,  and  in  particu¬ 
lar  for  the  higher  modes  (large  i/N  ratios)  where  N  >  9. 

3.  The  discussion  of  errors  centered  on  (l)  truncation  errors  from  the 
basic  central  difference  approximation  and  also  from  the  boundary- 
condition  representation  and  (2)  component  and  roundoff  errors.  Error 
considerations  as  well  as  the  ease  of  representing  boundaries  were  im¬ 
portant  in  choosing  to  approximate  the  first  spatial  derivative  rather 
than  the  second  or  higher  derivatives  and  also  in  the  use  of  staggered 
stations  in  the  calculation  of  successive  derivatives. 
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The  methods  proposed  in  this  report  are  relatively  insensitive 
to  component  or  roundoff  errors  for  N  <  1 6.  Larger  values  of  N  are 
not  recommended  because  of  the  general  computational  complexity  and 
also  the  increasing  size  of  component  or  roundoff  errors. 

h.  The  h4  and  symmetry  methods  are  nonphysical  in  the  sense  that  no  system 
composed  of  linear  bilateral  elements  can  be  found  that  is  represented 
exactly  by  the  equations  of  motion.  This  is  a  result  of  the  fact  that 
the  stiffness  matrices  are  not  symmetric  in  these  cases.  On  the  other 
hand,  the  h2  and  the  passive  circuit  methods  result  in  physically 
realizable  systems.  The  latter  two  methods  have  the  further  advantage 
that  the  resulting  modes  are  orthogonal,  i.e.,  a  transformation  to 
modal  coordinates  simultaneously  diagonalizes  the  mass  and  stiffness 
matrices. 


University  of  Michigan, 

Ann  Arbor,  Mich. ,  October  1961. 
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APPENDIX  A 

ALTERNATE  MATRIX  FORMULATION 


The  matrix  formulation  of  the  difference  equations  that  we  have  used  thus 
far  involves  the  solution  of  Eq.  (8l)  for  its  eigenvalues,  the  matrix  K  being 
obtained  after  three  matrix  multiplications.  Any  assumptions  with  regard  to 
boundary  conditions  or  constraints  are  put  into  the  individual  matrices  before 
the  matrix  multiplications  are  performed.  An  alternate  scheme,  particularly 
adapted  to  digital  computers,  enables  one  to  obtain  the  same  K  matrix  without 
losing  sight  of  the  subsidiary  assumptions. 

As  an  example  of  the  alternate  method,  consider  again  the  problem  of  the 
6-cell,  built-in  beam  using  symmetry  assumptions  at  the  ends.  Equations  (145) 
to  (148)  could  be  obtained  from  the  following  eight  equations: 

1  -27  27  -1  0 
0  1  -27  27  -1 
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These  equations  explicitly  separate  the  basic  difference  equations  and  the 
boundary  conditions. 
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APPENDIX  B 

TAPERED  CANTILEVER  BEAM 


The  analysis  that  has  been  made  in  the  main  body  of  the  report  has  been 
concerned,  with  the  lateral  vibrations  of  uniform  beams.  In  this  appendix  we 
will  indicate  the  sort  of  results  to  be  expected  when  higher-order  methods  are 
applied  to  nonuniform  beams. 


ANALYTICAL  SOLUTION 

The  case  of  a  beam  of  unit  width  and  with  uniform  mass  per  unit  volume  and 
a  linearly  decreasing  depth  has  been  analyzed  by  Siddall  and  Isakson. ^  Theoret 
ical  eigenvalues  were  obtained  for  the  case  where  the  depth  at  the  free  end  is 
one-fifth  that  at  the  built-in  end.  (See  Fig.  22.)  For  a  beam  of  unit  length, 
the  linear  density  is 

p  =  2(1  -  0.8x)  (Bl) 


and  the  bending  stiffness  is 

El  «  |(1  -  0.8x)3  (B2) 

The  cantilever  beam  is  built-in  at  x  =  0. 

Theoretical  frequencies  (rad/sec)  for  this  case  are 


Pi  =  2.47829 
e>2  =  9.08902 
p3  =  21.2953 


FINITE  DIFFERENCE  FORMULATION 

The  same  tapered  beam  was  studied  using  the  finite  difference  approach. 

In  a  manner  similar  to  that  for  the  uniform  beam,  the  natural  frequencies  were 
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calculated  using  various  boundary-condition  assumptions. 

For  a  nonuniform  beam,  nev  assumptions  muBt  be  made  concerning  the  method 
of  lumping  the  mass  and  stiffness.  Again  ve  assume  uniform  cell  size  ( h»l)  and 
end  the  beam  at  half- integer  stations.  The  total  mass  of  each  cell  is  lumped 
at  its  center.  Similarly,  the  total  effective  bending  stiffness  of  each  cell 
(i.e. ,  between  consecutive  half- stations)  is  lumped  at  its  center. 

Specifically,  the  mass  lumped  at  the  nth  station  is  (see  Fig.  22) 


“n 


(B3) 


Using  higher-order  differences,  we  obtain  the  acceleration  at  the  nth  station 
from 


yn 


V„.|  *  27V4  -  27Vn+i  ♦  W 


(Bit) 


The  lumped  bending  stiffness  is 


1  _  /  dx 

(El)n  *  J  El(x) 

n-1. 


(B5) 


resulting  in  the  bending  moment  equation 


«n 


Mn(e 

2kh 


n-fl 

2 


2T6n-i  +  279n+J.  "  en+^ 
2  2  2 


(B6) 


The  slope  and  shear  force  equations  remain  unaltered  from  the  case  of  a  uniform 
beam. 
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In  general,  the  boundary  conditions  are  imposed  in  the  seme  manner  as  for 
a  uniform  beam.  An  exception  occurs  in  the  passive  circuit  case  where  the  ef¬ 
fective  mass  m2  near  the  built-in  end  includes  the  reflected  value  (as  seen 
from  station  2)  of  the  actual  mass  mi.  A  similar  reflection  of  bending  stiff¬ 
ness  occurs  near  a  free  end. 

Also,  it  should  be  noted  that  the  constraint  equation 


yi  =  |ya  -  ^ya  (130) 

is  of  order  h4  for  the  case  of  a  nonuniform  beam  instead  of  order  h®  as  was  the 
case  for  the  uniform  beam. 


RESULTS 

Error  curves  for  the  tapered  beam  are  shown  in  Figs.  23-26.  Comparing 
these  results  with  those  for  the  uniform  beam,  one  finds  the  various  methods 
retain  many  of  their  characteristics  in  spite  of  a  general  reduction  of  accura¬ 
cy.  The  h2  method  shows  good  accuracy  for  the  fundamental  mode  while  the  higher- 
order  methods  (the  symmetry  and  passive  circuit  methods,  in  particular)  have 
better  accuracy  on  the  higher  modes.  The  h4  method  suffers  the  largest  loss 
of  accuracy  due  the  nonuniform  mass  and  stiffness  distributions. 

The  general  lowering  of  mode  frequencies  compared  to  the  uniform  beam  is 
probably  due  primarily  to  the  mass  lumping  method  which  results  in  a  slight 
positive  shift  in  the  center  of  mass  location  compared  to  the  continuous  beam. 
This  is  because  the  mass  is  lumped  at  the  center  of  the  cell  whereas  the  actual 
center  of  mass  for  the  cell  is  located  slightly  inboard  of  the  center,  i.e., 
toward  the  built-in  end.  Similarly,  the  compliance  (inverse  stiffness)  for  a 
cell  is  lumped  at  the  center  and  is  thus  somewhat  inboard  of  its  compliance 
centroid.  This  also  results  in  a  lowering  of  the  natural  frequencies. 
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TABLE  I 


ONE-SIDED  APPROXIMATIONS  OF  ORDER  h* 


A.  General  Case 


Sxl  A  24h 


77  [-  93yi  +  229y2  -  225y3  +  niy4  -  22y5]  +  p£\ 


64o  Sx5 


M  *  _Lr. 


[-  22yx  +  I7ya  +  9y3  -  5y4  +  ys]  -  -^h4 

1920  bx5 


SL '  se[- 22er  17ei +  +  ¥ 


71  h4  ase 
1920  bx5  ! 


B.  For  yA  =  0 
2 


S  7Trt3«T5yi  - 


1225y2  +  lt41y3  -  75y4l  +  -~1>4  |2 

128  ox  ^ 
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dx  U  840h  1920  bxs  2 


C.  For  yA  «  0,  9^-0 
2  2 
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TABLE  I  (Continued) 


D. 


For 

dx  A 


*  0 


r  *  919i  *  99Sa  -  99S  +  e2] 

dX  !  72n  2  2  2  2 


_L_h‘ 

1920 


dfe 

Sx=  ! 


76 


TABLE  II 


MODE  FREQUENCIES  OF  A  BUILT-IN  OR  FREE  UNIFORM  BEAM* 


Approximation  Method 

N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

3 

1 

22.373 

-  1.466* 
22.045 

3.831* 

23.230 

13.565* 

25.408 

4 

1 

22.373 

1.13511 

22.627 

1.850* 

22.787 

4.863 * 

23.461 

4 

2 

61.673 

-17.96111 

50.596 

-  1.516* 
60.738 

1.527* 

62.615 

5 

1 

22.373 

1.256* 

22.654 

0.796* 

22.551 

1.345* 

22.674 

2.700* 

22.977 

5 

2 

61.673 

-  9-357* 
55.902 

-II.O52* 

54.857 

0.028* 

61.690 

3.259* 

63.683 

5 

3 

120.90 

-27.833* 

87.246 

-24.607* 

91.155 

-  8.313* 
IIO.85 

-  6.907* 
112.55 

6 

1 

22.373 

1.059* 

22.610 

0.590* 

22.505 

1.095* 

22.618 

I.636* 

22.739 

6 

2 

61.673 

-  5.609* 
58.214 

-  3.755* 
59-357 

0.328* 

61.875 

2.679* 

63.325 

6 

3 

120.90 

-17.883* 

99.281 

-17.436* 

99.820 

-  3.201* 
117.03 

-  0.935* 
119.77 

6 

4 

199.86 

-34.099* 

131.71 

-23.131* 

153.63 

-14.000* 

171.88 

-13.274 * 
173.33 

*The  percentages  refer  to  the  error  in  frequency  as  compared  to  the  continuous 
beam. 
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TABLE  II  (Continued) 


Approximation  Method 

N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

7 

22.373 

0.858$ 

0.340$ 

0.912$ 

1.064$ 

1 

X 

22.565 

22.449 

22.577 

22.611 

7 

p 

61.673 

-  3.707$ 

-  0.751$ 

0.460$ 

2.006$ 

f 

59.387 

61.210 

61.957 

62.910 

7 

a 

120.90 

-12.333$ 

-  9.992$ 

-  1.472$ 

0.711$ 

P 

105.99 

108.82 

119.12 

121.76 

7 

u 

199.86 

-24.307$ 

-19.494$ 

-  6.915$ 

-  5*424$ 

1 

151.28 

160.90 

186.04 

189.02 

7 

c 

298.56 

-38.354 $ 

-23.617$ 

-18.412$ 

-18.000$ 

( 

P 

184.05 

228.05 

243.59 

244.82 

8 

T 

22.373 

0.697$ 

0.192$ 

0.764$ 

0.729$ 

22.529 

22.4l6 

22.544 

22.536 

8 

o 

61.673 

-  2.627$ 

0.091$ 

O.522 $ 

1.493$ 

60.053 

61.729 

61.995 

62.594 

8 

* 

120.90 

-  8.983$ 

-  4.806$ 

-  0.695$ 

1.092$ 

P 

110.04 

115.09 

120.06 

122.22 

8 

k 

199.86 

-18.063$ 

-14.140$ 

-  3.823$ 

-  2.086$ 

163.76 

171.60 

192.22 

195.69 

8 

C 

298.56 

-29.163$ 

-20.117$ 

-10.480 $ 

-  9.489$ 

P 

211.49 

238.50 

267.27 

270.23 

8 

6 

416.99 

-4l.4ll$ 

-24.933$ 

-21.823$ 

-21.569$ 

244.31 

313.02 

325.99 

327.05 
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TABLE  II  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

9 

1 

22.373 

0.572* 

22.501 

0.107* 

22.397 

0.653* 

22.519 

0.523* 

22.490 

9 

2 

61.673 

-  1-959* 

60.465 

0.243* 

61.823 

0.537* 

62.004 

1.129* 

62.369 

9 

3 

120.90 

-  6.832* 

112.64 

-  2.026* 
118.45 

-  0.281* 
120.56 

1.092* 

122.22 

9 

4 

199.86 

-13.905* 

172.07 

-  9.031* 
l8l.8l 

-  2.272* 
19502 

-  0.620* 
198.62 

9 

5 

298.56 

-22.773* 

230.57 

-16.342* 

249.77 

-  6.391* 

279.48 

-  5.078* 
283.40 

9 

6 

416.99 

-32.907* 

279.77 

-20.665* 

330.82 

-13.669* 

359-99 

-12.996* 

362.80 

9 

7 

555.17 

-43.707* 

312.52 

-26.471* 

408.21 

-24.504* 

419-13 

-24.337* 

420.06 

10 

1 

22.373 

0.474* 

22.479 

O.063* 

22.387 

0.550* 

22.496 

0.384* 

22.459 

10 

2 

61.673 

-  1.518* 
60.737 

0.217* 

61.807 

0.525* 

61.997 

0.866* 

62.207 

10 

3 

120.90 

-  5.368* 
114.41 

-  0.753* 
119.99 

-  0.033* 
120.86 

0.976* 

122.08 

10 

4 

199.86 

-11.018* 

177-84 

-  5.244* 
189.38 

-  1.401* 
197.06 

0.030* 

199.92 

10 

5 

298.56 

-18.214* 

244.18 

-12.165* 

262.24 

-  4.123* 
286.25 

-  2.726* 
290.42 

10 

6 

416.99 

-26.629* 

305.95 

-17.506* 

343.99 

-  8.947* 

379-68 

-  7-959* 
383.80 

10 

7 

555.17 

-35.865* 

356.06 

-21.478* 

435-93 

-16.449* 

463.85 

-15.972* 

466.50 

10 

8 

713-08 

-45.493* 

388.68 

-27.959* 

513.71 

-26.649* 

523.05 

-26.534* 

523.87 
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TABLE  II  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

ha 

h* 

Symmetry 

Passive 

Circuit 

11 

22.373 

0.402$ 

0.035$ 

0.478$ 

0.295$ 

1 

22.463 

22.381 

22.480 

22.439 

11 

o 

61.673 

-  1.213$ 

0.162$ 

0.499 $ 

0.678$ 

d 

60.925 

61.773 

61.981 

62.091 

11 

■z 

120.90 

-  4.334$ 

-  0.232$ 

0.116$ 

0.844$ 

J 

115.66 

120.62 

121.04 

121.92 

11 

J. 

199.86 

-  8.936$ 

-  2.867$ 

-  0.866$ 

0.314$ 

4 

182.00 

194.13 

198.13 

200.49 

11 

c 

298.56 

-14.871$ 

-  8.377$ 

-  2.767$ 

-  1.44o$ 

5 

254.16 

273.55 

290.30 

294.26 

11 

416.99 

-21.919$ 

-14.185$ 

-  6.098$ 

-  4.959$ 

o 

325.59 

357.84 

391.56 

396.31 

11 

f-f 

555-17 

-29.812$ 

-18.328$ 

-11.366$ 

-10.620$ 

7 

389.66 

453.42 

492.07 

496.21 

11 

8 

713-08 

-38.248$ 

-22.501$ 

-18.844$ 

-18.499$ 

440.34 

552.63 

478.71 

581.17 

11 

890.73 

-46.918$ 

-29.311$ 

-28.398$ 

-28.317$ 

9 

472.82 

629.65 

637.78 

638.50 
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TABLE  II  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

h2 

h* 

Symmetry 

Passive 

Circuit 

12 

22.373 

0.340* 

0.022* 

0.416* 

0.228* 

1 

22.449 

22.378 

22.466 

22.424 

12 

61.673 

-  0.992* 

0.114* 

0.467* 

0.538* 

d 

6l.06l 

61.743 

61.961 

62.005 

12 

■z 

120.90 

-  5.573 * 

-  0.025* 

0.199* 

0.720* 

II6.58 

120.87 

121.14 

121.77 

12 

j. 

199-86 

-  7.395* 

-  1.516* 

-  0.520* 

0.430* 

185.08 

196.83 

198.82 

200.72 

12 

c 

298.56 

-12.363* 

-  5-446* 

-  1.906* 

-  0.717* 

5 

261.65 

282.30 

292.87 

296.42 

12 

£ 

416.99 

-18.322 * 

-IO.856* 

-  4.295* 

-  3.132* 

D 

340.59 

371.72 

399.08 

403.93 

12 

r 7 

555.17 

-25.097 * 

-15.494* 

-  8.086* 

-  7.171* 

7 

415.84 

469.15 

510.28 

515.36 

12 

8 

713.08 

-32.468* 

-19.124* 

-13.590* 

-13.021* 

481.56 

576.71 

616.17 

620.23 

12 

890.75 

-40.206* 

-23.619* 

-20.910* 

-20.653* 

9 

532.60 

680.35 

704.48 

706.77 

12 

10 

1088.1 

-48.084* 

-30.508 * 

-29.850 * 

-29.786* 

564.93 

756.16 

763.32 

763.98 
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TABLE  II  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

h2 

h* 

Symmetry 

Passive 

Circuit 

15 

1 

22.375 

0.29511 

22.439 

0.013* 

22.376 

O.362* 

22.454 

0.179* 

22.413 

15 

2 

61.673 

-  0.827* 
61.163 

0.078* 

61.721 

0.431* 

61.959 

0.435* 

61.941 

15 

5 

120.90 

-  2.994* 
117.28 

0.041* 

120.95 

0.248* 

121.20 

0.6o4* 

121.63 

15 

4 

199.86 

-  6.219* 
187.43 

-  0.791* 
198.28 

-  O.285* 
199.29 

0.460* 

200.78 

15 

5 

298.56 

-10.430* 

267.42 

-  5-423* 
288.34 

-  1.550* 
294.59 

-  0.301* 
297.66 

15 

6 

416.99 

-15.526* 

552.25 

-  7.897* 
384.06 

-  5.105* 
404.05 

-  1.990* 
408.69 

15 

7 

555.17 

-21.379* 

436.48 

-12.656* 

484.91 

-  5-904* 
522.39 

-  4.914* 
527.89 

15 

8 

713.08 

-27.830* 

514.63 

-16.455* 

595.74 

-10.014* 

641.67 

-  9.278* 
646.92 

15 

9 

890.75 

-34.710* 

581.56 

-19.988* 

712.69 

-15.612* 

751.67 

-15.170* 

755.61 

15 

10 

1088.1 

-41.841* 

632.83 

-24.740* 

818.86 

-22.700* 

841.13 

-22.498 * 

843.26 

15 

11 

1505.5 

-49.054* 

665.02 

-51.564* 

895.52 

-31.075* 

899.71 

-31.027* 

900.30 

TABLE  II  (Continued) 


Approximation  Method 

N 

Mode 

Continuous 

h2 

h* 

Symmetry 

Passive 

Circuit 

14 

22.373 

0.255 $ 

0.009$ 

0.322 $ 

0.147$ 

22.430 

22.375 

22.445 

22.406 

14 

O 

61.673 

-  0.700$ 

0.052$ 

0.397$ 

0.355$ 

6l. 241 

61.705 

61.918 

61.892 

14 

x 

120.90 

-  2.548$ 

0.058$ 

0.28l$ 

0.521$ 

J 

117.82 

120.97 

121.24 

121.53 

14 

k 

199.86 

-  5.309 $ 

-  0.410$ 

-  0.130$ 

0.440$ 

189.25 

199.04 

199.60 

200.76 

14 

5 

298.56 

-  8.919$ 

-  2.124$ 

-  0.931$ 

-  0.060$ 

271.93 

292.22 

295.79 

298.38 

14 

c. 

416.99 

-13.314$ 

-  5.537$ 

-  2.28l$ 

-  1.266$ 

D 

361.47 

393.90 

407.48 

411.71 

14 

7 

555.17 

-18.407$ 

-  9.936$ 

-  4.408$ 

-  3-412$ 

i 

452.98 

500.01 

530.70 

536.23 

14 

8 

713-08 

-24.074$ 

-13-955$ 

-  7.531$ 

-  6.696$ 

541.41 

613.57 

659.38 

665.33 

14 

Q 

890.73 

-30.191$ 

-17.294$ 

-11.844$ 

-11.249$ 

y 

621.81 

736.69 

785.23 

790.53 

14 

10 

1088.1 

-36.623$ 

-20.908$ 

-17.434$ 

-17.085$ 

689.63 

860.61 

898.40 

902.19 

14 

11 

1305.3 

-43.231$ 

-25.837$ 

-24.263$ 

-24.107$ 

741.04 

968.06 

988.64 

990.63 

14 

12 

1542.1 

-49.867$ 

-32.482 $ 

-32.112$ 

-32.073 $ 

773.10 

1041.2 

1046.9 

1047.5 
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TABLE  II  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

hs 

h4 

Symmetry 

Passive 

Circuit 

15 

1 

22.373 

0.223 # 
22.423 

0.000# 

22.373 

0.282# 

22.436 

0.116#, 

22.399 

15 

2 

61.673 

-  0.602# 
61.302 

0. 034# 

61.694 

0.365# 

61.898 

0.293# 

61.854 

15 

3 

120.90 

-  2.200# 
118.24 

0.058# 

120.97 

0.289# 

121.25 

0.438# 

121.43 

15 

4 

199.86 

-  4.578# 
190. 71 

-  0.215# 
199.43 

-  0.020# 
199.82 

0.425# 

200.71 

15 

5 

298. 56 

-  7.714# 
275.53 

-  1.320# 
294.62 

-  0.646# 

296.63 

0.077# 
298. 79 

15 

6 

416.99 

-11.540# 

368.87 

-  3.801# 
401.14 

-  I.698# 
409.91 

-  0.794# 
413.68 

15 

7 

555-17 

-16.001# 

466.34 

-  7.536# 
513.33 

-  3.347# 
536.59 

-  2.392# 
541.89 

15 

8 

713.08 

-21.005# 

563.30 

-11.526# 

630.89 

-  5.765# 
671.97 

-  4.891# 
678.20 

15 

9 

890.73 

-26. 455# 
655.09 

-14.961# 

757.47 

-  9.128# 
809.42 

-  8.429# 

815.65 

15 

10 

1088.1 

-32.249# 

737.24 

-18.114# 

891.04 

-13.556# 

940.62 

-13.069# 

945.89 

15 

11 

1305.3 

-38.275# 

805.74 

-21.857# 

1020.0 

-19.084# 

1056.2 

-18.800# 

1059.9 

15 

12 

1542.1 

-44.413# 

857.23 

-26.859# 

1127.9 

-25.621# 

1147.0 

-25.498# 

1148.9 

15 

13 

1798.7 

-50.564# 

889.17 

-33.296# 

1199.8 

-33-007# 

1205.0 

-32.979# 

1205.5 

TABLE  II  (Concluded) 


Approximation  Method 

N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

16 

1 

22.373 

0.197# 

22.417 

0.000# 

22.373 

0.250# 

22.429 

0.098# 

22.395 

16 

2 

61.673 

-  0.524# 
61.350 

0.023# 

61.687 

0.33^# 

61.879 

0.245# 

61.824 

16 

3 

120.90 

-  1.911# 
118.59 

0.04l# 

120.95 

0.289# 

121.25 

0.380# 

121.36 

16 

4 

199-86 

-  3.993# 
191.88 

-  0.110# 
199.64 

0.050# 

199.96 

0.390# 

200.64 

16 

5 

298. 56 

-  6.736# 
278.45 

-  0.827# 
296.09 

-  0.439# 

297.25 

0.151# 

299.OI 

16 

6 

U16.99 

-10.096# 

374.89 

-  2.588# 
406.20 

-  1.273# 
411.68 

-  0.484# 
414.97 

16 

7 

555.17 

-14.030# 

477.28 

-  5.575# 
524.22 

-  2.576# 
540.87 

-  1.690# 
545.79 

16 

8 

713.08 

-18.471# 

581.37 

-  9.254# 
647.09 

-  4.483# 

68l.ll 

-  3.612# 
687.32 

16 

9 

890.73 

-23. 3*0* 
682.81 

-12.761# 

777.06 

-  7.141# 
827.12 

-  6.382# 
833.88 

16 

10 

1088.1 

-28. 562# 
777.32 

-15.817# 

915.96 

-10.671# 

971.99 

-10.082# 

978.39 

16 

11 

1305.3 

-34.046# 

860.86 

-18.946# 

1058.0 

-15.146# 

1107.6 

-14.748# 

1112.8 

16 

12 

1542.1 

-39.700# 

929.88 

-22.800# 

1190.5 

-20.556# 

1225.1 

-20.329# 

1228.6 

16 

13 

1798.7 

-45.438# 

981.40 

-27.809# 

1298.5 

-26.825# 

1316.2 

-26.725# 

1318.0 

16 

14 

2075.1 

-51.173# 

1013.2 

-34.018# 

1369.2 

-33.786# 

1374.0 

-33.762# 

1374.5 
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TABLE  III 

MODE  FREQUENCIES  OF  A  UNIFORM  CANTILEVER  BEAM 


Approximation  Method 

N 

Mode 

Continuous 

hs 

h4 

Symmetry 

Passive 

Circuit 

2 

1 

3.5160 

13.766# 

4.0000 

3.851# 

3.6514 

4.209# 

3.6640 

3 

1 

3.5160 

6.030# 

3.728 

1.812# 

3.5797 

1.257# 

3.5602 

3 

2 

22.034 

-  1.389# 
21.728 

3.980# 

22.911 

10.674# 

24.386 

4 

1 

3.5160 

3.356# 

3.634 

1.135# 

3.5559 

0.549# 

3.5353 

4 

2 

22.034 

1.035# 

22.262 

1.938# 

22.461 

5.328# 

23.208 

4 

3 

61. 697 

-17.941# 

50.628 

-  1.520# 
60.759 

1.483# 

62.612 

5 

1 

3.5160 

2.142# 

3.5913 

0.060# 

3.5181 

0.774# 

3.5432 

0.284# 

3.5260 

5 

2 

22.034 

1.153# 

22.288 

0.880# 

22.228 

1.402# 

22.343 

2.927# 

22.679 

5 

3 

61.697 

-  9.333# 
55.939 

-11.665# 

54.500 

0.023# 

61.711 

3.232# 

63.691 

5 

4 

120.90 

-27.841# 

87.242 

-23.956# 

91.937 

-  8.313# 
110.85 

-  6.898# 
112.56 

A 
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TABLE  III  (Continued) 


Approximation  Method 

N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

6 

1 

3.5160 

1.4818* 

3.5681 

0.028* 

3.5170 

0.555 * 
3.5355 

0.168* 

3.5219 

6 

2 

22.034 

0.976* 

22.249 

0. 613* 

22.169 

1.135* 

22.284 

1.770* 

22.424 

6 

3 

61.697 

-  5.589* 
58.249 

-  3.827* 
59.336 

0.324 * 
61.897 

2.663* 

63.340 

6 

It 

120. 90 

-17.885* 

99-277 

-17.443* 

99.811 

-  3.201* 
117.03 

-  0.935 * 
119.77 

6 

5 

199.86 

-34.099 * 
131.71 

-23.101* 

153.69 

-14.000* 

171.88 

-13.274* 

173-33 

7 

1 

3.5160 

1.086* 

3.5542 

0.017* 

3.5166 

0. 4l8* 
3.5307 

0.102* 

3.5196 

7 

2 

22.034 

0.794* 

22.209 

0.345* 

22.110 

0.944* 

22.242 

1.148* 

22.287 

7 

3 

61.697 

-  3.691* 
59.420 

-  0.772* 
61.221 

0.459* 

61.980 

1.997* 

62.929 

7 

4 

120. 90 

-12.333 * 
105.99 

-  9.992* 
108.82 

-  1.472* 
119.12 

0.711* 

121.76 

7 

5 

199.86 

-24.307* 

151.28 

-19. 479* 
160.93 

-  6.915* 
186.04 

-  5.424* 

189.02 

7 

6 

298. 56 

-38.354* 

184.05 

-23.620 * 
228.04 

-18. 412* 

243.59 

-18.000* 

244.82 
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TABLE  III  (Continued) 


• 

Approximation  Method 

N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

8 

1 

3.5160 

0.828* 

3. 51+51 

0.009* 

3.5163 

0.319* 

3.5272 

0.063* 

3.5182 

8 

2 

22.03L 

0.61+!+* 

22.176 

0.191* 

22.076 

0.790* 

22.208 

0.785* 

22.207 

8 

3 

61.697 

-  2.614* 
60.084 

0.086* 
61. 750 

0. 520* 
62.018 

1.488* 

62.615 

8 

4 

120. 90 

-  8.991* 
110.03 

-  4.797* 
115.10 

-  0.695* 
120.06 

1.092* 

122.22 

8 

5 

199.86 

-18.063* 

163.76 

-14.140* 

171.60 

-  3.823* 
192.22 

-  2.086* 
195.69 

8 

6 

298.56 

-29.163* 

211.49 

-20.117* 
238. 50 

-10.480* 

267.27 

-  9.489* 
270.23 

8 

7 

!+l6.99 

-4l.4ll* 

244.31 

-24.931* 

313.03 

-21.823* 

325.99 

-21.569* 

327.05 

9 

1 

3.5160 

0.648* 

3.5388 

0.006* 

3.5162 

O.236* 

3.5243 

0.043* 

3.5175 

9 

2 

22.03!+ 

0.531* 

22.151 

0.104* 

22.057 

0.667* 

22.181 

0.558* 

22.157 

9 

3 

61.697 

-  1.948 
60.495 

0.242* 

61.846 

0.534* 

62.028 

1.130* 

62.394 

9 

4 

120.90 

-  6.832* 

112.64 

-  2.026* 
118. 45 

-  0.281* 
120.56 

1.092* 

122.22 

9 

5 

199.86 

-15.905* 

172.07 

-  9.031* 
181.81 

-  2.272* 
195.32 

-  0.620* 
198.62 

9 

6 

298.56 

-22.773* 

230.57 

-16.342* 

249.77 

-  6.391* 
279.48 

-  5.074* 
283.41 

9 

7 

1+16.99 

-32.907 * 
279.77 

-20. 665* 
330.82 

-13.669* 

359.99 

-12.996* 

362.80 

9 

8 

555.17 

-43.707* 

312.52 

-26.471* 

408.21 

-24. 504* 
419-13 

-24.338* 

420.05 
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TABLE  III  (Continued) 


Approximation  Method 


N 

Mode 

Continuous 

h2 

h4 

Symmetry 

Passive 

Circuit 

10 

1 

3.5160 

0.509* 

3.5339 

0.003* 

3.5161 

0. 11+8* 
3.5212 

0.028* 

3.5170 

10 

2 

22.031* 

0. 1+1*0* 
22.131 

0.068* 
22.01*9  . 

0.567 * 
22.159 

0.409* 

22.124 

10 

3 

61.697 

-  1.509* 

60. 766 

0.217* 

61.831 

0.525* 

62.021 

0.864* 

62.230 

10 

120.90 

-  5.368* 
111*.  1*1 

-  0.753* 
119-99 

-  0.033* 
120.86 

0.976* 

122.08 

10 

5 

199.86 

-11.018* 

177.84 

-  5.241+* 
189.38 

-  1.1+01* 

197.06 

0.030* 

199.92 

10 

6 

298.56 

-18.214* 
21*1*.  18 

-12.165* 

262.21+ 

-  1+.123* 
286.25 

-  2.726* 
290.42 

10 

7 

1*16.99 

-26.627* 

305.96 

-17.506* 

31*3.99 

-  8.947* 
379.68 

-  7.959* 
383.80 

10 

8 

555.17 

-35.865 * 
356.06 

-21.1+78* 

1*35-93 

-16.447* 

463.86 

-15.972* 

466.50 

10 

9 

713.08 

-1*5.1*91* 

388.69 

-27.959* 

513.71 

-26.649* 

523.05 

-26.534* 

523.87 
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TABLE  III  (Concluded) 


Approximation  Method 


N 

Jlode 

Continuous 

u4  /1 - 1 - 

Passive 

n 

n 

symmetry 

Circuit 

11 

3.5160 

0.395* 

0.003 * 

0.054* 

0.023* 

1 

3.5299 

3.5161 

3.5179 

3.5168 

11 

22.03^ 

0.368* 

0.050* 

0.486* 

0.309* 

d 

22.115 

22.01*5 

22.141 

22.102 

11 

■j. 

61.697 

- 1.201+* 

0.162* 

0.499# 

0.676* 

3 

60.951* 

61.797 

62.005 

62.114 

11 

j. 

120.90 

-  l*.33l** 

-  0.232* 

0.116* 

0.835* 

4 

115.66 

120.62 

121.04 

121.91 

11 

c 

199.86 

-  8.936* 

-  2.867* 

-  0.866* 

0.315* 

182.00 

191*.  13 

198.13 

200.49 

11 

£ 

298.56 

-14.871* 

-  8.377* 

-  2.767* 

-  1.444* 

0 

25l*.  16 

273.55 

290.30 

294.25 

11 

H 

1*16.99 

-21.919* 

-14.185* 

-  6.096* 

-  4.959* 

( 

325.59 

357.84 

391.57 

396.31 

11 

8 

555.17 

-29.811* 

-18.328 * 

-11.366* 

-10.618* 

389.67 

453.42 

492.07 

496.22 

11 

713.08 

-38.250* 

-22.501* 

-18.844* 

-18.497* 

9 

1*1*0.33 

552.63 

578.71 

581.18 

11 

10 

890.73 

-1*6.918* 

-29.311* 

-28.398* 

-28.316* 

1*72.82 

629.65 

637.78 

638.51 

TABLE  IV 
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$ 

I 

°  o.  Mechanical  representation  of  the  simple  finite 

difference  approximation  to  a  uniform  beam. 


b.  Electrical  representation  of  the  simple  finite 
difference  approximation  to  a  uniform  beam. 


Fig.  3. 


o.  Even  symmetry  about  center  (symmetric  mode). 


i 


b.  Odd  symmetry  about  center  (antisymmetric  mode). 


Fig.  4. 
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R 


b.  Analog  computer  circuit  for  o  four-cell  cantilever 
beam  using  differences  of  order  h*  and  first- 
order  equations.  The  built-in  end  occurs  at 
station  1/2  and  the  free  end  at  station  9/2. 

Fig.  5. 


tr06-a 


t — WvJ 


b.  Odd  symmetry 


Fig.  6.  Anolog  computer  circuit  for  o  six-cell  beom  built-in  ot  both 
ends.  Differences  of  order  h*  ore  used.  Only  the  left  hoif  of  the 
beom  is  represented. 


a.  Even  symmetry 


b.  Odd  symmetry 

Fig.  7.  Analog  computer  circuit  for  a  six-cell  beom  thot  is  pinned  ot  both 
ends,  using  differences  of  order  h*.  In  (o.)  the  left  hoif  is  represented, 
while  in  (b.)  only  the  left  quorter  Is  represented,  due  to  symmetry. 
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using  the  central  difference  approximation  of  order 
h*  to  the  second  derivative  in  x. 


Fig.  10.  Symmetry  Assumptions  for  o  cantilever  beam.  Solid 
lines  show  assumptions  used  in  computations.  Dashed  lines 
indicate  deviations  giving  a  physically  consistent  set  of 
assumptions. 
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Fig.  II.  Passive  circuit  onalogy  for  fht  basic  h4  approximation. 


Fig.  12.  Passive  circuit  ana 
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for  a  built-in  end  at  station  1/2. 


ERROR  IN  PERCENT 
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Fig.  15.  Built-in,  h4  approximation. 
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Fig.  16.  Cantilever,  h4  approximation. 
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Fig.  17.  Built-in,  using  symmetry. 
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Fig.  19.  Built-in,  passive  circuit. 
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Fig.  21.  Pinned,  h4  approximotion. 
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